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Abstract 



We derive explicit formulae for a set of constraints for the Einstein 
equations on a null hypersurface, in arbitrary dimensions. We solve 
these constraints and show that they provide necessary and sufficient 
conditions so that a spacetime solution of the Cauchy problem on a 
characteristic cone for the hyperbolic system of the reduced Einstein 
equations in wave- map gauge also satisfies the full Einstein equations. 
We prove a geometric uniqueness theorem for this Cauchy problem in 
the vacuum case. 
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1 Introduction 

An Einsteinian spacetime is a pair (V, g) with V a manifold and g a Lorentzian 
metric solution on V of Einstein's equations. These are a geometric system 
of quasilinear second order partial differential equations which express, in 
vacuum, the vanishing of the Ricci tensor of g. Two isometric spacetimes 
are considered as the same. A fundamental problem is the determination of 
generic solutions from initial data, i.e. the Cauchy problem. It is a badly 
posed problem from the analyst point of view in arbitrary local coordinates 
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on V since the characteristic determinant of the PDE system is identically 
zero. It is well known that if initial data are given on a spacelike manifold 
M the problem splits into constraints to be satisfied by the geometric initial 
data, i.e. the induced metric and second fundamental form of M as subman- 
ifold of V, and an evolution system. The latter must be well posed and must 
exhibit propagation of the gravitational field with the speed of light, that 
is, must be a hyperbolic system with characteristic cones the null cones of 
the looked-for spacetime metric. The simplest way to obtain such a system 
is to suppose that the coordinates satisfy so-called harmonicity conditions, 
or, more generally, to introduce a preassigned metric g, called target metric, 
which permits to write the Ricci tensor as the sum of two tensorial opera- 
tors, one of which is a hyperbolic operator acting on g, called the reduced 
Ricci tensor, and the other a homogeneous first order differential operator 
acting on a vector H, called wave-map gauge vector, which vanishes when 
the identity map is a wave map from {V,g) onto {V,g). When M is space- 
like, classical theorems of analysis show existence and uniqueness of solutions 
of so-reduced Einstein equations and, due to the Bianchi identities, of the 
original Einstein equations if the initial data satisfy the constraints. The 
case where the initial manifold is null has analogies with the spacelike case 
but also important differences: First, the induced metric is degenerate, and 
unconstrained in the regions where r (as defined below) has no zeroes. Next, 
the second fundamental forms defined on a spacelike and on a null manifold, 
for which the normal is also tangent, have very different properties. Finally, 
null initial data on a light cone, or on two-intersecting null hypersurfaces, 
determine the solution in one time direction only, past or future. 

A complete understanding of this problem is still lacking, even in space- 
time dimension four. The most exhaustive studies are for the case of two 
intersecting null surfaces [SISlIISlinillSllSSlElllSlllHlEn]; compare, in different 
settings, [11 [21 [32]. The most complete construction of equations satisfied 
by initial data has been given by Damour and Schmidt [17], and the most 
satisfactory treatment of the local evolution by Rendall [16]. The problem 
with data on a characteristic cone presents new mathematical difficulties due 
to its singularity at the vertex, and only partial results have been obtained 
before in [8l|20l|2ll[2ll[17] . 

The object of this work is to present a treatment of the Einstein equations 
with data on a characteristic cone in all dimensions n + 1 > 3. We proceed 
as follows: 

Though the equations are geometric and the final results coordinate in- 
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dependent, it is useful to introduce adapted coordinates to carry-out the 
analysis. We take a C°° manifold V diffeomorphic to R""*"^, and we con- 
sider a cone Co in V with vertex O & V and equation, in coordinates 
compatible with the C°° structure of V 

y' = r, r:={Y^ {y^V^ . 

i=l,...,n 

We consider the Cauchy problem with data on Co for the Einstein equations 
with unknown a Lorentzian metric g, assuming that Co will be a charac- 
teristic cone of the metric g and the lines = r, y = d, where the c* are 
constants, its null rays. It is well knowi]0that the characteristic cone of a C^'^ 
Lorentzian metric admits always such a representation in a neighbourhood 
of its vertex. We review in section 3 an existence theorem which applies to 
the reduced Einstein equations in wave-map gauge with Minkowski target 
reading in these coordinates 

g = -(dy'f+ (dyr- (1-1) 

i=l,...,n 

We introduce in section H] what we call adapted null coordinates, singular on 
the line r = 0, in particular at the vertex O of Co, but C°° elsewhere, by 
setting 

x° := r — y^, x'^ := r, 

and defining x^, A = 2, ...n, to be local coordinates on the sphere S*""^. In 
coordinates x" the trace 'g on Co of the metric g we are looking for has the 
form 

g = gQQ{dx'^y + 2i'Qdx^dx^ + 21/ ^dx^dx^ + 'gj^^dx^dx^ . 

Remark that the question, whether or not x^ is an affine parameter on the 
null rays x"^ = c"^, depends on derivatives tranversal to Co of the spacetime 
metric g, which are usually not considered as part of the initial data for 
characteristic Cauchy problems. 

The adapted null, but singular at the vertex, coordinates are used to 
solve "wave-map gauge constraints" satisfied by 'g. 

^It is known, using normal coordinates, that a null cone in a C^'^, Lorentzian spacetime 
always admits such a representation in a neighbourhood of its vertex; the null geodesies 
which generate it are represented by the lines where r only varies. This is guaranteed to 
hold only in a neighbourhood of the vertex, as there can be caustics. 
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In section we review the standard argument, that the Bianchi identi- 
ties imply that if g satisfies the reduced Einstein equations with source a 
divergence-free stress energy tensor, then the vector H satisfies a homoge- 
neous hyperbohc system; it vanishes in the future of Co if its trace H vanishes 
on Co ■ 

We show in sections [H] to [TT] that if = if the initial data 'g satisfies a 
set of + 1 equations which we call the wave-map gauge constraints. These 
constraints read as a hierarchical system of ordinary differential equations 
along the light rays, singular at the vertex O, if one uses the adapted null 
coordinates x". We write this complete system for a general g and gener- 
alized wave gauge, in arbitrary dimensions n + 1 > 3. We integrate them 
successively under natural limit conditions on the unknowns at O. We study 
briefiy in section 16.21 the case when the degenerate metric g induced on Co 
(i.e. the dependent quadratic gAB) is prescribed. 

In section [7]in order to have an evolutionary equation for r we set, as many 
authors before us, = ^2^7^^, with 7 an arbitrarily given dependent 
metric on S"""^. The first wave- map gauge constraint can be written in a form 
which involves the two unknowns and f2. Its general solution is obtained 
by the introduction of an arbitrary function k. We study in particular the 
case K = which leads to the Raychaudhuri equation for r for which we prove 
global existence for a small \a\ which depends only on the given 7. A simple 
integration determines then f2, hence and we are back to the equations 
for i/q, VAi'goQ with given g. We remark that the equation for z/q (for k = 0) 
implies that the vector I is parallelly transported along the null ray by the 
connection of a spacetime metric in wave-map gauge satisfying the Einstein 
equations. In sections [H [9l [10] and [11] we establish, and integrate, the other 
constraints determining va and ^qo- ^ theorem in section [T2] summarizes our 
analysis of the wave-gauge constraint equations. A uniqueness theorem is 
proved in Section [13] 

A major question left open by our work is the description of the largest 
class of unconstrained initial data which lead to solutions of the wave-map 
gauge constraints such that the components in coordinates of the trace 'g 
satisfy the (non trivial) initial conditions given in section [311] for the existence 
theorem for quasilinear wave equations. The problem is that the wave-gauge 
constraint equations determine the components of 'g in the coordinates, 
and these components are linked with the components in the coordinates 
by linear relations which are singular at the vertex. We simply note here 
that initial data which are Minkowskian in a neighbourhood of the vertex 



6 



are easily seen to be in the class where the existence theorem holds; see 
also [12] for more general family of data. We plan to return to this problem 
in a near future. 

2 Definitions 

2.1 Ricci tensor and harmonicity functions 

The Ricci tensor of any pseudo Riemannian metric is given in local coordi- 
nates by, 

A A A A ^ 

Rap ■= dxT^i^ - daTf^^ + r^^^r^^ - r^^r^;^, dx ■■= -g^, (2.1) 

with the Christoffel symbols 

:= 9^''[f^, [/W, := ^{^agp^, + d/^gaf, - d^Qap). (2.2) 

The Ricci tensor satisfies the identity 

Rap = R% + \{9axdpT^ + gpxdaT''), (2.3) 

where Ricc^^\g), the reduced Ricci tensor, is a quasi-linear, quasi- diagonal 
operator on the components of 

rI^^ ^ -^g^f^dxd.g^p + f[g, dgU , (2-4) 

and f[g, dg]^^ is a quadratic form in the first derivatives dg of g with coeffi- 
cients polynomial in g and its contravariant associate. 
The r^'s, called harmonicity functions, are defined as 

r'^ := g^^Vl^ . (2.5) 

The condition = expresses that the coordinate function satisfies the 
wave equation in the metric g. 
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2.2 Wave-map gauge 

The harmonicity functions are coordinate dependent and only defined locally 
in general, whether in space, or time, or both. The wave-map gauge, which 
we are about to define, provides conditions which are tensorial. A metric g 
on a manifold V will be said to be in g-wave-map gauge if the identity map 
V ^ V is a. harmonic diffeomorphism from the spacetime (V, g) onto the 
pseudo-Riemannian manifold {V,g), with g a given metric on V. Recall that 
a mapping / : (V, g) — )■ (V, g) is a harmonic map if it satisfies the equation, 
in abstract index notation, 

□r := 9^'{di^r - nAf^ + d^rd^rv-j = o . (2.6) 

In a subset in which / is the identity map defined by /"(a;) = the 
above equation reduces to H = 0, where the wave-gauge vector H is given 
in arbitrary coordinates by the formula 

■= c/"^r^^ - , with := ^"^f , (2.7) 

where f are the Christoffel symbols of the target metric g. See [S] for a 
more complete discussion of the concepts and results in this section. 

The following identity has been proved to hold, with D the Riemannian 
covariant derivative in the metric g O page 163], 

Raf^ = Rif + ^{gaxD^H'' + g^xDaH""), (2.8) 

where R^^^^ (g), called the reduced Ricci tensor of the metric g in ^-wave- 
map gauge, is a quasi-linear, quasi-diagonal operator on g, tensor-valued, 
depending on g: 

Rif = -Ig^'DxD.gap + f[g, DgU , (2.9) 

where f[g,Dg]ai3, independent of the second derivatives of 5^, is a tensor 
quadratic in Dg with coefficients depending upon g and g, of the form (see 
formula (7.7) in chapter VI of |9j) 

P!:;''^''{g)Dpg^,sDagx>. + Ig^nsapRx'^f. + gppRx'a,} , (2.10) 
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with R the Riemann curvature tensor of the covariant derivative D. We 
will frequently restrict ourselves to the case in which the target metric is 
the Minkowski metric 1], and then denote by D the covariant derivative. In 
this case, and if using coordinates such that the Minkowski metric takes the 
canonical form (II. ip . the reduced Ricci tensor in wave-map gauge coincides 
with the one in harmonic coordinates. 

We emphasise that, unless explicitly stated, our computations are valid 
for a general g. 

Our main results below assume that W takes the form (12. 7p . However, 
several results apply to a large class of W^s of the form|§ 

■.= g''^t^^p + W\ (2.11) 

where is a vector which may depend upon g, g and possibly some other 
fields, but not upon the derivatives of g; the relevant restrictions are pointed 
out in (HSD-dZZD, (K^ - (K^ and (^MM-<MM- The reduced Ricci tensor 
becomes then 

RiT'' = Ri? + \ {gc.xDpW'^ + gp^Kw>^) . (2.12) 

However, unless explicitly indicated otherwise we assume that W is identi- 
cally zero. 

Another interesting generalization (see, e.g., [H] and references therein) 
has been inspired by numerical simulations: if one uses the decomposition 
(12. Sp . the identities = are only obeyed to some finite precision and 
shows a generic tendency to deviate from zero. Attempts to cure that 
have been made by introducing constraint damping terms [30], changing the 

( H) 

choice of the reduced Ricci tensor -R^^ to 

Rff)^ + \(^{n^9l3\ + n^g^x - -^^^nxgap)H^ , (2.13) 
or, equivalently, the reduced Einstein tensor 5*^^ to 

+ ^^^agpx + npga,x + (p - l)nxgap)H^ , (2.14) 

^H.Friedrich [25] introduced generalized harmonic coordinates by adding arbitrary func- 
tions to the harmonicity conditions. 
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where n'^ is a vector field and e is a small positive constant which controls 
the rate of damping of the gauge conditions. (As shown in [50] the constant 
p must also be positive to have damping.) We will show that the damping 
terms are consistent with our analysis. For definiteness we will assume that 
has been prescribed, though certain more general situations can easily be 
incorporated into our scheme. 



3 Characteristic Cauchy problem 

The Einstein equations in wave-map gauge with source a given stress-energy 
tensor T, 

^fp = -i^a^^DxD^go,^ + f[g, Dg]^,^ = p^p , p^^ := T^/j - -^^^dap , 

form a quasi-diagonal, hyperquasi-lineaiG system of wave equations for the 
Lorentzian metric g. The Cauchy problem for such systems with data on a 
spacelike n-manifold Mq is well understood, the Cauchy data are the values 
of the unknown on Mq and their first transversal derivatives. When Mq is 
not spacelike in the spacetime {V^g) which we are going to construct, the 
problem is more delicate. It is known since Leray's work (see [37]), @ that 
the Cauchy problem for a linear hyperbolic system on a given globally hy- 
perbolic spacetime is well posed if Mq is "compact towards the past" ; that 
is, is intersected along a compact set by the past of any compact subset of 
V . However the data depend on the nature of Mq and the formulation of 
a theorem requires more care. In the case where Mq is a null hypersurface, 
except at some singular subsets (intersection in the case of two null hyper- 
surfaces, vertex in the case of a null cone) the data is only the function, not 
its transversal derivative, with some hypotheses which need to be made as 
one approaches the singular set. 

In this article we concentrate on the case of the light cone, though most 
of the calculations of our equations apply to any null hypersurface. 

''That is, the principal second order terms are diagonal and their coefficients depend 
on the unknowns but not on their derivatives 

^See [551fi^ for a treatment of generalized solutions of a linear wave equation with data 
on a achronal Lipschitz section of a spacetime with compact spacelike sections. 
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3.1 The Cagnac-Dossa theorem 

To prove the local existence of solutions of Einstein equations with data on a 
characteristic cone we use a wave-map gauge, and an existence theorem for 
solutions of quasi linear wave equations with such data. 

The proof of an existence theorem for such a characteristic quasilin- 
ear Cauchy problem is inspired by the Leray's idea of the linear case, ap- 
plied to the characteristic cone and linear wave equations in Cagnac [6] (cf. 
also Cagnac |5] and Friedlander [22]); extended to the quasilinear case by 
Cagnac [7]. The most complete results appear in Dossa's thesis, the second 
part of which is published in abbreviated form in [12]. One considers quasi- 
diagonal, quasi-linear second order system for a set v of scalar functions w^, 
/ = 1, . . . , iV, on R"+i of the form 

A^>^{y, v)dl^v + /(y, v,dv) = Q, y= {y^) G R"+\ A, /i = 0, 1, . . . , n > 2 , (3.2) 
v = {v^). dv = {^) dlv = {^), f = {f), / = !,..., AT. (3.3) 

The initial data 

V := v\co = (f) (3.4) 

are given on a subset, including its vertex O, of a characteristic cone Cq- 
Throughout this work a bar over an object denotes the restriction of that 
object to Co- 

Cagnac and Dossa assume that there is a domain U C R"^^ where Co is 
represented by the following cone in R""*"^ (compare Figure B?T] below) 

Co:={x'^r-y' = 0}, := (v')'- 

1=1,. ..,n 

The initial data is assumed to be defined on the domain 

C^ := Co n {0 < t := 2/° < T} . (3.5) 

They denote 

Yo := {t := ?/° > r} , the interior of Co , := Yo D {0 < < T} . 

(3.6) 

They also set 

:= Co n {y^ = t} , diffeomorphic to 5""^ , (3.7) 
Sr := Yo n {y^ = t} , diffeomorphic to the ball 5""^ . (3.8) 



11 



We will use the following theorem given in the first part of Dossa's thesis: 
it assumes some more differentiability of the data than the theorem in |T2], 
but it is simpler to apply to the Einstein equations whose initial data must 
satisfy wave-map-gauge constraints, and is sufficient for us here. 

Remark that these results assume more regularity from the data on the 
cone than the regularity obtained for the solution, a constant fact in chara- 
teristic Cauchy problem already seen in other contexts. 

Theorem 3.1 Consider the problem i\3.2\) - l!i3.4\ )- Suppose that: 

1. There is an open set U d R"^"^ x R^, Yq C U where the functions 
A^^ are (7^™+^ in y and v. The function f is C^™ in y eU and v E W and 
in dv G R("+^)^. 

2. For {y,v) E U x W the quadratic form A^^^ has Lorentzian signature; 
it takes the Minkowskian values for y = and v = 0. 

3. a. The function cf) takes its values in W . The cone Cq is null for the 
metric A^^{y, cj)) and (f){0) = 0. 

b. (f) is the trace on Cq of a (7^™+^ function in U. 

Then there is a number < Tq <T < +oo, Tq = T if (p is small enough in 
(j2m+2 jiQj-Yyi^ gy^(,]i fij^Qf fiiQ problem 113. ijj)-([5^ has one and only one solution 



V m Yq" , such that 



1. If m > ^ + 1, V e K^'+^Y^") n in particular \dv\ is 
bounded. 

2. If m = oc, V can be extended by continuity to a C°° function defined 
on a neighbourhood of the origin in R^+^. 

The spaces K"^{Yq) and F"^{Yq) are Banach spaces of sets of functions 
on Yq which together with their time and space derivatives of order less or 
equal to m admit a square integrable restriction to each St and for which, 
respectively, the following norms are finite 



1=1,. ..,N -'^ 0<|fc|<m 



M\f^{y,^)-= E sup t 2 V ||5V||l2(5,) . 

1=1,. ..,N 0<|fc|<m 

The Euclidean metric, e := J^ii^y^)^, is used to define the measure on St 
and as usual k denotes a multi- index, k := {ko, ki, . . . , kn), the derivation 
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ko + ki + . . . + kn- 



:= {do)'°idi)''^ . . . (dn)'" , with 9„ := — . (3.9) 

3.2 Einstein equations in the wave-map gauge 

We know that the wave-map gauge reduced Einstein equations on a manifold 
V are tensorial equations under coordinate changes, so that any coordinates 
can be used. Note that the principal part of the wave-map reduced Einstein 
equations is independent of the target manifold, and so the Einstein equations 
on R""^^ in wave-map gauge are of the form (13.21) for an unknown h, when 
we set g = rj + h and work in the y coordinates where the Minkowski metric 
takes the standard form 

^ ^ -(dy'f + J^Wf- 

i 

As an application of Theorem 13. II we obtain (see also [20j in space-dimension 
n = 3): 

Theorem 3.2 (Existence for the wave-gauge reduced Einstein equations.) 
Let 'g = T] + h be a quadratic form on Cq such that the components h^y 
in the coordinates y^ satisfy the hypotheses of the Existence Theorem \3.1\ 
Then, if the source p is of class C^™ in Yq , there exists Tq > such that 
the wave-gauge reduced Einstein equation^ R^^^ = p^f^ admit one and only 
one solution on Yq° , a Lorentzian metric g'^^'' = rj + h, with h satisfying the 
conclusions of that theorem. 

The reader should keep in mind that a solution of the Einstein equations 
in wave-map gauge as above will be a solution of the original Einstein equa- 
tions if and only if the wave-map gauge vector H defined by the constructed 
spacetime metric is such that Rice = Ricc^^\ condition satisfied in particular 
iiH = 0. 

The following theorem is a straightforward adaptation of a theorem proved 
long ago by one of us [21] for spacelike Cauchy data. 

^We use abstract index notation when it helps formulate properties of geometric objects. 
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Theorem 3.3 Let g^^^ be a Lorentzian metric, solution on Yq of the 
Einstein equations in wave-map gauge Sj^^ = Tap- Then g^^^ is a solution 
on Yq of the full Einstein equations Sap = Tap if the wave-gauge vector 
vanishes on Cq and the source T satisfies the conservation law V aT"'^ = 0. 

Proof. The identity (12 .81) imphes (indices raised with g) 

gap ^ gaP(H) ^ 1 + _ g^f^D^H^). (3.10) 

Hence the equations in wave-map gauge Sj^'' = Tap and the Bianchi identities 
imply that H satisfy the quasidiagonal linear homogeneous system of second 
order equations 

Vab^'H^ + Vab^H'' -V^baH"" = , (3.11) 

whose principal terms are wave equations in the metric g since Vq, 
V^DaH°' is at most first order in H . If g is (7^, H is C^, and an energy 
inequality applied to this linear equation implies easily that = in Yq if 
H ■= H\rT = Q. ■ 

When the support of the initial data is a spacelike manifold Mq the van- 
ishing of H is guaranteed when the constraint equations {Sap—Tap)n^\Mo = 
are satisfied by the initial data, where is the field of unit normals to Mq in 
the space-time one seeks to construct. One of the main goals of this work is 
to present a method to construct initial data on the light-cone which ensures 
the vanishing of H. 

4 Null hypersurfaces, adapted coordinates 

The obtention, and solution, of equations to be satisfied by initial data to 
ensure the vanishing of H is simpler in coordinates adapted to the geometry 
of the null initial manifold. 

4.1 Adapted coordinates 

Let Mq be a hypersurface in R"+^ which will be a null submanifold of the 
spacetime (y,g) with V some domain of R""*"^. Mq is generated by geodesic 
null curves, called rays. In a manner classical for null surfaces we choose 
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coordinates x°' so that Mq is given by the equation x° = 0, and on Mq the 
coordinate X IS cl parameter along the rays, denoting by i the tangent vector 
We assume that the subspaces Sj-i : {a;^=constant, a;°=0} are spacehke 
and diffeomorphic to the same n — 1 manifold E, except possibly for Sq which 
reduces to a point in the case of a characteristic cone. We denote by x"^ local 
coordinates on E. We have = 0, = 1, = 0. 

The assumption that Mq, x'^ = 0, is a null surface for g is equivalent to 
saying that 'Qh = 0. The covariant vector n := grada;°, with x^ = the 
equation of Mq, is a null vector normal and tangent to Mq with components 
no — 1, rii — ua — 0. By uniqueness of null directions tangent to a light cone 
we have also £a — and hence, using that i"' — S^, g^^ — 0. Then the trace 
on Mo of the spacetime metric reduces in the coordinates to (we put an 
overbar to denote restriction to Mq of spacetime quantities) 

g :— g\xO=o = ^oo(^^°)^ + 2vQdx^dx^ + 2vAdx^dx^ + g^^^dx^^dx^ , (4.1) 

where 

'^o-=9qi, iyA-=9QA, (4-2) 

We remark that the are the non zero components of the quadratic form 
g induced by g on Mq by the identity map. They define an -dependent 
Riemannian metric on E 

h ■= gAsdx^dx^, A,B = 2,...,n. (4.3) 

The following identities hold on Mq, because g^^ and g""^ are inverse 
matrices. 

r^f^^O, u':=f'^-, (4.4) 

:tAB _ ~AB ~AB 



We denote 
then 



g = g , with g the inverse matrix of g^^^. (4.5) 

:= g^^UA (4.6) 
g^' ^ r ^ -i^ygoo + {A'^^^A. (4.7) 



Remark 4.1 We consider a null adapted coordinates condition for the trace 
Ij but, in contradistinction with Bondi-Sachs and other authors, we do not 
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assume these coordinates for the spacetime metric. Such an assumption is a 
gauge condition which can be considered as the analogous for the character- 
istic Cauchy problem to the zero shift and constant lapse of the n + 1 classical 
decomposition, gauge which is not well adapted to the well posedness of the 
evolution problem. 

In Appendix |A] we collect formulae useful for explicit calculations, such 
as the trace on Mq of the Christoffel symbols of g, etc. 

See [2H1I221I311I1I] for various useful results concerning null surfaces. 

4.2 Characteristic cones 
4.2.1 General properties 

It is no geometric restrictioij^ to assume that in a neighbourhood of its vertex 
the characteristic con^ of the spacetime we are looking for is represented in 
some admissible coordinates y := (y") = y*, i = 1, . . . , n) of R"'^^ by the 
equation of a Minkowskian cone with vertex O, 

r-y' = 0, r:={J2iyyy- (4-8) 

i 

Given the coordinates we can define coordinates on R"+i adapted to 
the null cone Co as we did for a general null surface by setting 

a;0 = r - = r, x^ = /i^(^) , (4.9) 

r 

with x"^ local coordinate^ on S*""^. The null geodesies issued from O have 
equation x° = 0, x"^ =constant, so that is tangent to those geodesies. On 
Co (but not outside of it in general) the spacetime metric g that we are going 
to construct takes the form fl4.ll) . that is, such that g-^^^ = and ^^^^ = 0. 

We emphasize that our assumption 'g is given by (14. ip is no geometric 
restriction on a metric g to have such a trace, however ^qq, z/q, i^a are not 
invariant under an isometry of spacetime with leaves Co invariant, they are 
gauge-dependent quantities (see sections 1^31 and [T3l) . 

^See footnote [T] and details in section 23] 

cone is a topological manifold but it is not difFerentiable at its vertex 
^ They can be angular coordinates, see e.g. [131 Chapter V, Section 4], or stereographic 
coordinates, as in Christodoulou |14j . 
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We compute the relation between the components of a tensor T in the 
coordinates y and x using fl4.9p and its inverse: 

= x^ - y' = x^e^(x^), with ^(e^^ = 1 . (4.10) 

i 

If the components of a spacetime symmetric tensor T in the coordinates 
are denoted Tq,/3 and if in the coordinates y" they are denoted Tq/j, then 

the identities for the change of coordinates of tensors, Tx^ = Ta^||x||^, give 
the identities 

■ 96* 
^00 = Tqoi Tqi = —Too ^ TQiQ\ Tqa = ~T^'q^ ' (4-11) 

Tn = Too + 27o,e' + T^Q'Q\ T^a = lotr^ + TjfQ^^ , (4-12) 

Conversely, 7^ = f^f^T^/? gives 

Too = Too, Zk=-(^oo+Toi)e*-ToA^, (4.14) 

T,, = (Too + 2Toi + Tn)e^e^ + (To^ + Ti^)(e*^ + e^-^) + T^B^7- , 
— ay-' oy^ oy^ oy^ 

(4.15) 

with 

dx^ -1 /I 

where the /if 's are C°° functions of the x^ on any subset of S*""^ where the 
x'^'s are admissible local coordinates. 
One checks, using the identities 

i i 

with 

SAsdx'^dx^ = Sn-i , the metric of S"'~^ , 
that, when T = t], 
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4.2.2 Case g = t], the Minkowski metric 



It is natural to take as given metric g the metric of a model spacetime such 
as Minkowski, or de Sitter, or anti-de Sitter. While most our formulae will 
be completely general, the analysis will mainly be concerned with the case 
where the metric ^ is a Minkowski metric rj given by the formula written 
above in the introduced coordinates and in the adapted null coordinates 
x". The Riemannian curvature of the Minkowski metric t] is zero. The 
non zero Christoffel symbols of rj are in our coordinates x", with S^q the 
Christoffel symbols of the metric s, 

r?A - , f^c - sic , (4.16) 

= -x^sab , = -x^sab ■ (4.17) 

Equalities and identities assuming given metric g = i] and = 5'"'^f will 

be denoted with symbols = and =, respectively when ambiguous. 
We have 

W'^-x'g^'^SAB^w', (4.18) 

i 2t'^Tfc + r^'ric = -4^°^^ + g'^'s^c . (4.19) 

4.2.3 Limits at the vertex 

We set g = T] + h. The condition hapiO) = of Theorem 13.21 can always 
be satisfied by choice of an orthonormal frame for the natural frame of the 
coordinates y°' at the vertex. Since the coordinates x are singular for = 
the behaviour near x^ = of the components hap in x coordinates is obtained 
only by considering limits. As explained above, we can, and will, choose 
coordinates on Co such that hn = 0, i.e. Cq : x" = is a null cone for 
g, and hiA = i.e. the vector ^ := is on Cq a null vector. Then the 
components of h are 

Qx"^ 

hm = hoo, Ki = -{hQQ+hm)Q'-hoA-Q— with /iqi := uq-I, Ka ■= t^a , 

= (hoo + 2hoi)&e^ + 7^oA(e^-- + e^--) + Hab 



dy^ dy^ dy"^ dy- 
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We see that the condition hap{0) = is equivalent to the following 
conditions in the coordinates x°': 

lim(l + g^o) = lim(z/o - 1) = lim(r-VA) = lim {r-^TiAB) = . (4.20) 

) — ^0 r— >-0 1 — >0 r— >-0 

4.2.4 A lemma 

For further use we note the following observation: 

Lemma 4.2 // a spacetime function f is such that on Cq in the coordi- 
nates it holds that 



lim dif = 



then it also holds 



lim dof — . 

r=x^—yO 



Proof. One has the trivial identity 

If / is in a neighbourhood of O, dif tends to a limit, a number at O, 
hence the above equation implies 



limSoZ + aiO' = , 

r->-0 



condition which can be satisfied for all if and only if Oj = and limr_».o dof — 
0. Therefore 

lim dof = — hm dof = . 

r— >0 r->0 



4.3 The afRne-parameterization condition 

The vector field ^ '■— tangent to the null rays in Mq, obeys the geodesic 
property 

PV^F = , (4.21) 

with 

rn = rfi = and T\ = i^\d^vo-]f^i) . (4.22) 
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If we impose the condition = 0, then the vector i is parallelly transported, 
and is said to be an affine parameter on the rays; this condition gives an 
equation involving uq, a metric coefficient which will appear in our first wave- 
map gauge constraint. However, we stress that the equation F^^^ = involves 
also a derivative transversal to Mq, and thus cannot be made to hold just by 
a gauge choice of the coordinate on the cone. We will see later how we 
can circumvent this problem in the wave-map gauge. 

4.4 Null extrinsic curvature 
4.4.1 General properties 

Let Mq be a null hypersurface with a field of null tangents i. The null 
extrinsic curvature at x G Mq is defined (see, e.g., [28]) as the bilinear form 
with components Vq,£/3 acting on the quotient of the tangent space to Mq at x 
by the direction defined by i, i.e. equivalence classes of tangent vectors of the 
form X =X+cC. with T^Mq, c an arbitrary number. Indeed, the action 
of the bilinear form on a pair of such tangent vectors, Va^/jX^y^, depends 
only on the equivalence class, that is in our coordinates on the components 
and Y^, hence it is defined by the components Xab •= Va^b of the 
bilinear form. Using = 5" and la '■= Qap^i = ^o^a have 

Xab = -^AB^o = \9i'gAB- (4-23) 

We denote by 

Xa^ := ^^^Xac = rf, (4.24) 

the mixed, x^-dependent, 2-tensor on S^^^ deduced from the null second 
fundamental form. We define its trace 

r := Xa^ = g^^'xAB = dii\og ^detgj,), (4.25) 
and its traceless part 

:= —^A^^ and we set \a\^ := (Ta^ctb^ . (4.26) 

n — 1 

See [28||29] for an analysis of the null second fundamental form through 
the Weingarten map. 
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4.4.2 Harmonicity functions 

In adapted coordinates, and using tlie notation above, the harmonicity func- 
tions = g^^T'^^ reduce on Mq to 

r° ^ ^^-r;, ^ 2z.°rS, + ^^^r^^ ^ A'^'d^i - r) , (4.27) 

^^'^T^^ = 2z.°r^ 2^^T^, +J^Ti^ (4.28) 
- 2u%\b^ + f^, (4.29) 

= ^^^ri^ = ^^t;, + 2z.t;, - 2u'u^T\, + r^'r^^ (4.30) 

= (^°)'5i^oo + ^''^°(^^ + 5i^o - ruo) + 2(z/°) V^(-9ii.A + ly^AB) 
+ A^'^Vbz/a - (4.31) 

= - d,t' + t'l^'ild^i - dii^o - Tuo) + A^'^Vbz/a - V^^. 

(4.32) 

We have defined 

T^:=r^tic, (4.33) 
with T^fj being the Christoffel symbols of the metric g^s- We shall also use 

Ti := g,^r = Uof, (4.34) 
:= gAB^^ ^ '9a,J^ , (4-35) 
and similarly for components of W and H with subindices. 

4.4.3 Vertex limits 

We set 

9ab = r\sAB + Jab) ■ (4-36) 

We have seen in Section 14.2.31 that it is no geometric restriction for smooth 
metrics to assume 

lim {r'^gAB - Sab) = , i.e. lim/^^ = , 

r— !>0 r— >0 
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Then 



limrV^'^s^'' ■ (4.37) 

r— >-0 



Recalling that Xa" = ^"ai using the relation between connections in 
different frames gives 



V ^ = f - ^-^ , 



Using 



Qya q^a Qx'i-^ dy"^ dx'^ dx^ 



dy^ ^ dy^ y^ ^ dx'-^ d dy°' 1 dx^ dy"^ 1 

dx^ ' dx^ r ' dy"^ dx^ dx^ r dy'^ dx^ r ^ ' 

we find 

Therefore if the coefficients r*Q and V''^^ are bounded for < r < a, the 

same property holds for for Xa^ ~ ^^A' ^^^^ ~ ^'^'^ "^A*" 

Xa*" ~ ^^'^a''"- These quantities are then also continuous on each null ray. 
However the limits when r tends to zero are in general angle dependent. 

We have already said no geometric restriction for a given spacetime metric 
to choose for it normal geodesic coordinates in a neighbourhood of a point, 
here O. Then the Christoffel symbols vanish at O and if continuous are such 
that 

dx^ d a?/" _ 1 dx^ dy'^ _ 1 



r->o dy'^ dx^ dx^ r dy"' dx^ r ^ 



Hence 



c 

- mil u / 

See further results in next section. 



lim ^ = lim a a = 0. (4.38) 



4.5 Boundary conditions in coordinates normal at the 
vertex 

In the following sections we will give explicit expressions for the wave-map 
gauge constraints. To study their solutions we will need to know the be- 
haviour of the unknowns at the tip of the light cone, aiming at finding so- 
lutions of the constraints which satisfy the Cagnac-Dossa hypotheses. The 
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purpose of this section is to describe this behaviour in coordinate systems 
useful for the problem at hand. The analysis here is also useful for proving 
geometric uniqueness of solutions. 

Consider a smooth space-time (V, g). Let O E V and let Co be the future 
light-cone emanating from O. Let T be any unit timelike vector at O, and 
normalize null vectors £ at O by requiring that g{i,T) = —1. The parallel 
transport of i defines an affine parameter, denoted by s, on the future null 
geodesies s i— 7f(s) with 7^(0) = O and with initial tangent £. Let {z'^), 
/i = 0, . . . , n be a normal coordinate system centred at O with T = d^o, see 
e.g. [Mllin], or [9l Chapter 12, Section 7]. In those coordinates the future 
light-cone emanating from O is given by the equation 

n 

Co = {z° = \z\} , where z := {z\...,z'') , |zp := J^^z')^ . 

i=l 

As is well known, in normal coordinates at O, 2; = 0, the Christoffel symbols 
vanish at O. Hence for a C^'^ metric we have daQfiuiO) = 0, and so, for small 

\z\ := -|- \z\, 

\9i.v - Vl + \4dagt.A < C^kP , (4.39) 

for some constant C. 

In the coordinate system (x^) = = m, = r, x"^), A = 2, . . . , n, where 

u=\z\- z^ , r =\z\ , (4.40) 

and where the x'^'s are any local coordinates on S*""^ parameterizing the unit 
vector z/\z\, the trace of the metric g on Co takes the desired form (14. ip as 
long as the metric, and the light-cone, are smooth; assuming smoothness of 
g, this will always be the case in some neighbourhood of the tip O. 
Equation (I4.40p shows that 

z^ 

dz^ = dr — du , dz^ = — dr + dAZ^dx"^ , 

r 

which allows us to translate the estimates (14.390 to the asymptotic behaviour 
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of the objects of interest near r = 0: From 
g^Jz^'dz" = {7]^^ + 0{\z\^))dz^'dz'' 

= (-1 + 0{\z\^)){du - dry + 0{\z\^){du - dr) ( -dr + dAz'dx' 



r 



+ {6{ + 0{\z\^)) (^jdr + dAz'dx-^^ (^^dr + dsz^dx^^ 

(-1 + 0{\z\^)){duY + (2 + 0{\z\^))dudr + 0{r\z\^)dudx^ 
+0{\z\^){dry + 0{r\z\^)dx^dr + r^ {sab + 0{\z\^)) dx^dx"" 



we obtain, at m = 0, for small r, 



= -1 + 0(r2) , drVoo = 0{r) , d^= 0(r) , dAQoo = 0{r^) , (4.41) 

uo = l + 0{r^) , druo = 0(r) , d^gp i = 0(r) , dA^o = 0{r^) , (4.42) 

UA = 0(r3) , drUA = 0(r2) , dj^ = 0(r2) , dAUo = 0(r3) , (4.43) 

9AB = rHsAB + 0(r2)) , drigAB - r'sAB) = 0(r3) , (4.44) 

duQAB = 0(r3) , dAjgA B - r'sAB) = 0{r^) , (4.45) 

gn=0 , drgu = , dugu = 0{r) , dAgu = , (4.46) 

9iA = 0, drg,A = , d^^g^ = 0{r^) , dAgiA = • (4.47) 

One also has associated second-derivatives estimates, 



dudr gAB = O(r^) , d'^.{9AB - r'^SAB) = O(r^) , dAdr{gAB - r'^SAB) = 0{r^) , 

(4.48) 

etc. From fl4.44p and f l4.48p we obtain 

Xa"" = Ma + 0{r) , hence r = ^ + 0(r) , = 0(r) , (4.49) 
as well as 9,(r - ^) = 0(1) , Oat = 0{r) , (4.50) 
draA^" = 0(1) , dcaA"" = 0(r) . (4.51) 

Note that f l4.49p - fl4.51l) will hold in any coordinate system which coincides 
with the normal coordinates z^ on the light-cone. This is due to the fact that 
the vectors dr and Sa are tangent to the light-cone, which implies that the 
quadratic form g^. = gAB^x^dx^ is intrinsically defined on the light-cone, 
independently of how the coordinates are extended away from the light-cone, 
and from the fact that the matrix g^^ in flA.2p is the inverse of 'gAB- 
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7^2 




Figure 4.1: The cross-section of the hght-cone Cq] Cq is the shaded blue 
region. Two generators and 7^^ are also shown. 

4.6 The light-cone theorem 

A result closely related to our analysis here is the light- cone theorem, proved 
in [11], which reads as follows: Let s be an affine parameter as defined at 
the beginning of Section H75l Let denote the (n — 1) -dimensional surface 
reached by these geodesies after affine time s: 

= {7,(s)} C Co , (4.52) 

where the vectors i run over all null future vectors at O normalized as above; 
see Figure 14. 1[ Note that in this subsection 14.61 and only here, we shall 
slice the null cone Co and its interior using an affine parameter s, and not 
coordinate time as in the rest of the article. 

We denote by Cq the subset of the light-cone covered by all the geodesies 
up to affine time t: 

C'o = Uo<.<tS, . (4.53) 

Note that 7^(s) might not be defined for all s. Further, some of the S^'s 
might not be smooth. However, there exists a maximal Sq > such that 
is defined and smooth for all < s < Sq. Our considerations only apply to 
that last region. 
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It is proved in [TT] that, assuming the Einstein equations with a cosmo- 
logical constant and with sources satisfying the dominant energy conditions, 
the areas of the S^'s are less than or equal to the corresponding areas in 
Minkowski, de Sitter, or anti-de Sitter space-time. Furthermore, if equality 
holds at some S2, then on Cq we have 

— n — 1 

0"AB = = Til , T = . 

r 

(This situation will be referred to as that of the Null-Cone Theorem (NCT).) 
It is further shown in ^Jj that, under suitably stronger energy conditions, 
equality implies that the metric is that of the model space on the domain of 
dependence of Cq . The proofs of those facts provide a non-trivial illustration 
of the formalism developed here, as specialized to the simpler problem treated 
in [n]. 

5 Constraints and gauge preservation 

The obvious analogue on a null submanifold Mq of the spacelike constraints 
operator is Sapl^ ■, where ^ denotes the field of null normals to Mq normalized 
in some arbitrary way. Derivatives of the metric in SnH^^ transversal to 
the light-cone appear only at first order. Some of them □ cancel magically 
between the various terms contributing to Sa/si^, and those that remain 
can be expressed in terms of H. So, in the explicit form of Sa/si^, one 
can replace every occurrence of d^goi, OqQqa and 'g^^d^gAB by H^, Wa, and 
terms containing only derivatives along Mq. We define n+1 operators Ca{H), 
a = 0, . . . , ra, by adding all the terms involving H in Sapd-^ ■ One can then 
define n + 1 operators Ca by whatever remains; thus the C^'s coincide with 
Sap^^ when vanishes. Explicit formulae for Ca are given in fl6.13p . f l8.24p 
and in ffTOTj) below, while can be found in f lCTj) . and f lTaBH]) . 

We will prove the following theorem, which is the key element of our 
analysis of the Cauchy problem for the Einstein equations on a light-cone: 

Theorem 5.1 1. The operator Sapd-^ on a null submanifold Mq can he 
written as a sum, 

^Compare [17l|46] in space-dimension three. 
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where vanishes when H = 0, while the operator Ca depends only on 
the values 'g on Mq of the spacetime metric, on the choice of the null 
vector i, and on W , which depends on the chosen target space of the 
wave-map gauge. The operators Ca will be called Einstein-wave-map 
gauge constraint operators. 

2. In adapted null coordinates: 

a) The operators Ca lead to a hierarchy of ordinary differential oper- 
ators for the coefficients ofg along the generators, all linear when 
the first constraint Sai3^°''^^ = Ta/ji'^i^ has been solved. 

b) The operators Ca together with the wave-gauge reduced Einstein 
equations lead to a hierarchy of homogeneous first order ordinary 
linear differential operators along the generators for the compo- 
nents Ha if the spacetime metric g satisfies on Mq the reduced 
Einstein equations. 

Theorem 15.11 will be proved in Sections (MHl and [101 

A consequence of Theorems 13.31 and 15.11 is the following: 

Theorem 5.2 A Lorentzian metric g^^\ solution of the Einstein equa- 
tions in wave-map gauge S'^J^ = Ta/3 in Yq with VaT"^ = 0, and taking an 
initial value 'g on Cq, is a solution of the full Einstein equations Saj3 = Tap 
if and only if'g satisfies the constraints Ca = Tapi^ . 

Proof. Theorem 15.11 gives the following identities, with Ca a linear and 
homogeneous first order differential operator along the null vector i for the 
vector H, 

Sapi^ = Sl^pf + \{gaxDpH^ + gpxDjf' - gapDjf'y^ = Ca+ Ca. (5.1) 

The "only if part of the theorem results inmediately from the identity 
(15. ip when the metric g is a solution of the full Einstein equations and is in 
wave gauge, since then only Ca remains in that identity. 

The "if part will be proved later by showing that H = is the only 
solution, for metrics which are uniformly near the tip of the cone, of the 
equations 

^{gax'^^ + ^/3aA^ - 9ap'KH~^)i^ = Ca (5.2) 
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which result from the identity (15. ip when Ca = T^id^^ and 5^^^ = Tap- ■ 

The question of local geometric uniqueness of solutions is addressed in 
Section [T31 



6 The first constraint 

6.1 Computation of 'Rn = Sn = Sa/s^''^^ 
The component Rn can be separated as 

D _ p(l) , d(2) 

where is linear in first derivatives of the Christoffel symbols and R^i^ is 
quadratic in them. They are given by, after a trivial simplification, 

nf^ ^ 9or?, + dATf, - d,r% - d,rf^ , (e.i) 
i?S? ^ r?,(r; + Tl, + r^j + tI,{t% + t\, + vf^) + Tf,{T% + r^, + r^^) 

pO pO 9pl pO pl pl 9pl P"'^ 9pO P^ P^ P^ ^'fi 9"^ 

10-*- 10 ^-"^ 10-*^ 11 ^ 11-*^ 11 ^-"^ lA-"^ 11 lA'^ 10 IB-*- lA- l^-^J 

We must take care when taking derivatives transversal to the cone, i.e. 
do, that our coordinates are valid only on the cone. We will then use the 
trivial identity 



dx^^ = r'{dx[fi, /37]) + (5A^7"^)[/i, H- (6.3) 

mly is 

9iA = 0, 



In R^^^ only T^j^ is differentiated transversally to Ca- We have, since g^^ 



doT^n = ^'^'d.dogu + {d^uo - -dogu)dog^', with dog^' = -{uydogn, 



By using also 



(6.4) 

9ir?o^-9i(iz/°^). (6.5) 



rfi = 0, Tt = lg^^d,gAs=:T (6.6) 
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and the harmonicity function fl4.27p . we get 



= {iy')''^dogn do9n + ^diu'^dogn - d^r (6.7) 
= ^(fi + rf - ^z/°9ii/o(ri + r) - d,r. (6.8) 

(2) 

The part R^^ depends only on the values of the Christoffel symbols on 
Co- Using T^-^ = Til = T^^ = and trivial simplifications we find that 

^(2) ^ pi /pO , N _ pO pO _ pA =5 

^11 — 111-*- 10 "T IAJ lO-*- 10 IB^ lA ■ 

(2) 

In the chosen coordinates, R^^ reduces to 

R^^^ ^ -l(ri+r)2 + iz/°9iZ.o(ri+r)+z/°9iZ/or-^r(ri+r)-XA''xB^ • (6-9) 
Adding (16. 8 p and (16. 9p we obtain 

Rn = -diT + u'd,uoT-^r(r, + T)~XA''xB^ (6-10) 
^ -^^r + T\,T-XA''xB^, (6.11) 

6.2 The Ci constraint operator 

By definition of the wave-gauge vector H we have Fi = Wi + Hi, and hence 
(I6.10p decomposes as 

Rn = Ci + Ci, (6.12) 

with 

Ci := -diT + f z/°9iZ/o -kwi + t)] t - \af - , (6.13) 

\ 2 / n — 1 

where we have separated Xa^ trace- free and pure trace parts (see (I4.26P ). 
and 

A := -i^iT . (6.14) 

As announced (16.130 involves only the values of the metric coefficients 
on the light-cone; equivalently, no derivatives of the metric transverse to the 
light-cone occur there: 

Wi = 2Vli + z/o^^^f% ^ -v^x^g^^'sAB , (6.15) 
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where we have assumed that the target metric takes the adapted form in the 
same coordinate system, so that F'j'i = ^'^^ — 0- The Einstein equation 
Rii = Til in wave-gauge provides, in this sense, a constraint equation Ci = 
Til for the metric components g^^. 

The constraint equation Ci = Tu contains as unknowns only the compo- 
nents ^jiQ and z/q if it is so of Tu. A simple strategy is then to prescribe 
(compare Bondi et al. [31j) and use the definition (I4.23p to compute Xab^ 
hence also and r. The first constraint reads then as a differential first or- 
der equation for i^q, linear if Tu is independent of z/q since Wi is linear in uq. 
(Recall that we are assuming = unless explicitly indicated otherwise.) 
The solution will lead to a Lorentzian metric as long as uq is positive. 

However, the equation is singular wherever r vanishes, as the resulting 
ODE for z/q involves inverse powers of r. For this reason it is of interest to 
look for alternatives, where r is computed from the constraint, rather than 
prescribed in advance. Following |46] we will prescribe only the conformal 
class of g. The wave- map gauge constraint deduced from f l6.13p is then an 
equation for uq and the conformal factor Q^. We can prescribe arbitrarily 
z/q and then determine Q. We can also, generalizing an idea of Damour 
and Schmidt, impose to uq to satisfy a well chosen differential equation con- 
taining an arbitrarily given function k. We treat in detail the case k = 0, 
which implies that for the obtained solution uq the vector i will be parallelly 
transported, in other words r will be an affine parameter, in the resulting 
space-time. 

7 Solution of the Ci constraint for given a 

The operator Ci relates the three functions r (which, via equation fl4.25p . 
essentially describes the evolution of the volume element of the sections E), 
z/q and |crp := aA^crs^- We recall the following well known fact: 

Lemma 7.1 The tensor a is determined by the conformal class of the in- 
duced quadratic form g . 

Proof. To see this, let us write 
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with 7 a degenerate quadratic form on Cq such that = = 0. Then 

Xab = ^^^^l7AiJ + lAB^diQ , 

and thus 

X^^ ^ \l''^d,^^s + 5^ d, log n , (7.1) 
hence the trace-free part of Xa^ is 

^^a"^ = -7 ^i7ab - ^5;^^ . (7-2) 

where 7^] denotes the positive definite dependent quadratic form on E 
with components 7^5. We see that the traceless tensor cr is independent of 
the conformal factor, hence depends only on the conformal class of ^. In 
particular a vanishes if g is conformal to a quadratic form independent of 
r := x^. ■ 

If 7 and its first derivatives satisfy the vertex limits spelled out for g in 
section H75l then limrlcrl = 0; we say that a degenerate quadratic form 7 on 

r-J-O 

Co, with 7^]^ = 7j^^ = 0, is admissible if it is on Cq — O, i.e. for r > 0, 
and such that |(Tp is for r > 0, hence bounded for finite r > 0. Given a 
the constraint Cu = Tu appears as a relation between the functions r and 
uq. Since it involves radial derivatives of both r and (which can actually 
be grouped as (9i(i/°r)) we could prescribe one of them and integrate for the 
remaining field, or else provide an additional differential equation for either 
of T or z/q and integrate simultaneously the coupled system of the constraint 
and this new equation. In the remainder of this section we show how to solve 
the constraint by prescribing z/q, either explicitly (section 17. ip or through a 
differential condition (rest of section [7]). 

7.1 Prescribed uq 

Suppose the function uq is arbitrarily prescribed, then the constraint equation 
becomes a differential equation for r. It is convenient to introduce the scalar 
function (recall that denotes the restriction of g to E) 

VdetSn-i/ \detSn-iJ 
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so that 

r = (n — 1) 9i log v?, or diip = (p. (7.4) 

n — 1 

The normahzation of p has been chosen to have (p = r for Minkowski. Using 
this variable the constraint reads 

- "^nV? + u^'diuo - -VTi di<p = — <p , 7.5 

\ z I ip J n — \ 

to be integrated outwards with initial data V'(O) = and d\p){p) = 1. As 
already mentioned, W \ contains ^p nonlinearly, and in principle Tn could 
also depend on p). In general this scheme could be considered for a larger 
class of gauge 

^i = W^i(7AB,¥',^o,r,x^), (7.6) 

and sources of the form 

T\\ = Tn (source data, 7^6. ^iTAB. di(p, z/q, diUo, r, x^) , (7.7) 

where di denotes derivatives tangential to the light-cone, and by "source 
data" we mean non gravitational data, for example fields determined from 
characteristic initial data for scalar, or Maxwell, fields. The wave-map gauge 
condition f l4.18p is clearly of the form (17. 6p . In Section [72] we show that both 
scalar and Maxwell fields lead to a coefficient Tn of the energy-momentum 
tensor compatible with (17. 7p . 

7.2 Differential equation for uq 

The choice made by Rendall is to assume that is an affine parameter 
along the null rays, in other words that the vector i = is parallelly trans- 
ported along the null rays by the connection of the spacetime he constructs, 
i.e. T^^ = 0; equivalently di^o = (Fi + r)i/°/2. Now, this last equation 
contains a derivative transversal to the light-cone, which is not part of the 
characteristic initial data. Extending to the cone an idea of Damour and 
Schmidt [17] concerning two intersecting surfaces, in anticipation of the fact 
that our solution will satisfy i^i = 0, we could impose to z/q to satisfy the 
equation 

diiyo = liWi + T)iyo, (7.8) 
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which imphes, modulo Hi = 0, that t\i = 0. When uq satisfies (I7.8P the 
constraint Ci = Tn reduces to a Raychaudhuri type equation for the only 
unknown r 

diT+ ^ + |a|2 + Tn = 0. (7.9) 

n — 1 

More generally all solutions of f l7.9p can be obtained by introducing an arbi- 
trary function k, and solving the pair of equations 

diT-KT + ^ + |a|2 + Tn = 0, (7.10) 

n — 1 

whose only unknown is r when and Tn are known, and 

diUo=^(Wi + T)uo + KUo . (7.11) 

When uq satisfies this equation and Hi = 0, then = n. 

Once r is determined we can use (17. 4p to obtain and finally (17. lip to 
compute z/q. 

Remark 7.2 The equations derived here would be dramatically simplified 
if one simultaneously imposes z/q = 1 and k = 0; see, e.g., [16]. However, 
these two conditions together with the wave-gauge condition, which is the 
cornerstone of our analysis, would impose undesirable geometrical restrictions 
on the initial data. 

Equation (I7.1QP is, if Tn does not depend on r or a Riccati differential 
equation on each null ray and hence can be rewritten, precisely using the 
variable y?, as a linear second-order equation 

- 9iiV2 + k9iv? = -^^^-^ t^V^ , (7-12) 
n — 1 

to be integrated outwards with initial values ^{O) and diipiO). We have as- 
sumed chosen an admissible 7, hence |(Tp continuous for r > 0. We assume 
the same holds for Tn. Remark that for a continuous tensor T, i.e. with con- 
tinuous components T^i, in the coordinates y", we will have Tn continuous, 

but lim,f._^o ^11 a function of angles in general since it holds that 

limTn = roo(O) + 2T,o(0)^ + T,,(0)^ . 
r-*-o — r — - 
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When |(Tp + Tii is continuous for r > standard ODE theory guarantees 
that a solution with given initial values exists globally. However, a positive 
definite metric is only obtained from the positive part of the solution. 
The relevant initial conditions are ^p{0) = and diip{0) = 1, so ip is initially 
tangent to ip = r. 

We consider the = 0, that is, is an affine parameter. Assuming 

Til > 0, the equation satisfied by p shows that it is a concave function of 

on each null ray wherever positive, and hence there are two possibilities: 
either p> is a monotone increasing function for all real r, with < (p < r 
and < r < (n — l)/r, or else there is a first local maximum, at which 
di(p = and hence the expansion r also vanishes there. This is related to 
the formation of outer-trapped surfaces on the cone Co- Once a maximum 
has been reached, ip will necessarily vanish for some larger value of r. 

We now turn to a direct analysis in terms of r, which allows stating results 
in a more geometric way. 



7.3 Solution of the Raychaudhuri equation 

We continue to use a Minkowski target and we make the choice k = 0, so 
that x^ will be an affine parameter along the null rays. Equation (I7.10p then 
reads as a Raychauduri equation 

diT + ^—T^ + \a\^ + Tii = . (7.13) 

72—1 

This is a first order ODE for r when |crp := crA^<^B^ and Tn are known. 



7.3.1 NCT case 

When |crp + Tu = 0, the equation admits the solution corresponding to the 
Minkowskian conj^: 

n — 1 , , 

r, = ^. (7.14) 

The value |(t| = further imposes 

1 2 
Xa^ = ^^A, i-e- dig^B = -^9ab ■ (7-15) 
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It is the only solution such that r ^ tends to zero at the vertex of the cone. 
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With our choice of frame of coordinates at the vertex, the solution is the 
Minkowskian solution 

9ab = {x^fsAB , (7.16) 
as used in the null-cone theorem of [TTl. 



7.3.2 General case, a global existence theorem 

We denote by r and di by a prime. In all that follows we write, and 
solve, differential equations in r, with constant initial values (mostly zero) 
for r = 0. We do not write explicitly the dependence on the other coordinates 
x^, though it occurs in the solutions and in the coefficients. 

1. In a neighbourhood of r = we define a new function y by 

n — 1 
y ■= • 

T 

Equation (17. 13p becomes 

y' = l + ^fV, f:=\a\' + Tu. (7.17) 
n — 1 

In agreement with Section r4.4.3[ we seek a solution such that y{0) = 0. The 
equation implies that y is increasing and y > r. 

We assume as before that ^^y^^ continuous and bounded by a number 
A^. Then y exists, is of class C^, is unique, and is bounded by the solution 
of the problem 

z' = 1 + A^z^, z(0) = , 
as long as that solution exists. The solution is 

z = A-^ tBxi{Ar) . (7.18) 

Hence z is defined, and bounded, as well as all its derivatives, for < 
r < a, for any a < A~^^. 

For < r < a < A"'^, z is such that 



r < z <r + A^r^ 

We have defined ip as 



11 — \ 

ijj ■=r^- r, Tr, = , (7.19) 
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and hence we have, since r < y < z, 

1,11111 1 A'^r .2 z^, X 

< ^ < < = T— ^ < A'r . 7.20 

n — I r y r z r r + A-^r-^ 1 + ^^r^ 

That is 

< < (n- 1)AV . (7.21) 
2. For large r we use the decay of p. Using the definition fl7.19p we 



obtain 



^' + = + f . (7.22) 

r n — 1 



u = -— + r^f > , u:=r^ij . (7.23) 



This gives 

1 

n — 1 r'^ 

Hence u is an increasing function of r, and both u and ijj are positive as 
u{0) = 0. 

For r > a we replace the problem to solve by the integral equation 

uir) = ^ [ + r p'f{p)dp + , (7.24) 

with Ua := u{a) = a^ip^a). By fl7.18p and the inequality y < ^ for r < a we 
have 

,,, n — 1 .x/l A 

ip{r) < Trj{r) — = (n - 1) - 



z{r) \r tan(y4r) 

hence 

We assume that r^p is integrable for r G [a, oo), and we set 



oo 



J a 

The solution u of the integral equation f l7.24p exists and is bounded by a 
solution V of the equation 

= / ^^dp + Ca , 



n-l Ja P 
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as long as such a solution v exists; equivalently, as long as the differential 
equation 

1 



v' 



n — Ir"^ 

admits a solution v with v{a) = Ca- The general solution of the above 
equation is 

11 in — l)r 
- = — h c , I.e. V = — . 

V [n — l)r 1 + (n — l)rc 

It takes the value Ca for r = a if and only if 

'^^'^ Ca, i.e. c = Ca := - - — (7.26) 



1 + {n — l)ac Ca {n — l)a 

The function v remains positive and bounded if 

l + {n- l)rca > , (7.27) 

hence v is defined and bounded for all r if Ca > 0, i.e. Ca < {n — l)a; that 
is, when 

Ua + Ba<in- l)a . (7.28) 
It follows from f l7.25p that this last inequality will hold when 

Ba<in-1)- 



tan(y4a) 



In the case where = for < r < a it holds also that u{r) = in that 
interval we have Ua = and Ca = Ba- Condition f l7.28p reduces to 

Ba<{n- l)a . (7.29) 

3. Assume that r > exists in the interval r G (0,6] and denote Tb '■= 
T{b). If for r > b, p = 0, the equation for r reduces, for r > 6, to 

t' + ^-r^ = , (7.30) 
n — 1 



with initial value 



Ti — 1 

r(6) = n, < Tfe < ro(6) = — — . (7.31) 



37 



The solution is 

1 r 1 



(7.32) 



T n — 1 Tb n — 1 
Therefore 

^ (n - l)Tb ^ (n-l) 

(n-l) + (r-6)n r + 4' ^ ' 

with 

4 := (n - l)r^^ - 6 > 0. (7.34) 

Hence, for large r, 

r = I^(l-^ + ...). (7.35) 
r r 

We have proved the following theorem. 
Theorem 7.3 The equation for r deduced from the first constraint, 

r' + ^r^ + f = 0, 
n — 1 

with := |(T^| + Til continuous and r^f^ integrable in r for r e [0, oo), has 
a global solution T(r) > 0, and the function 

n — 1 

ip := r 

r 

is of class if: 

1. We assume that there exists a e (0, oo) such that it holds 

TT 1 

A<—, with := sup -/^ . (7.36) 

2a 0<r<a n — 1 

In the interval <r < a it then follows that 

n — 1 , , . 

T > z{Ar) , 

r 

with 

z(x) := < 1. (7.37) 

^ ^ tanx - ^ ^ 

2. In the interval a < r < oo we assume that 

poo 

/ r^ f^ dr < {n - l)a z{Aa) . (7.38) 

J a 
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In this interval we then have 

n — 1 



T > 



r + ka 
with 



poo 

ka = in-l){Ta-a-''Ba)-'-a, Ba= r^fdr. (7.39) 

J a 

3. Regardless of point 1., if a = = Tu for r > b, and if := T{b) > 0, 
then the solution exists for all r > b and it holds that 

"""^ k['^ -.= (71-1^7' -b>0. 



Remark 7.4 If := |crp + Tu has compact support {a < r < b} with 
a > 0, it follows from (I7.28p that f l7.38p can be replaced by 



f r^f^dr <{n- l)a 

J a 



which will be satisfied if, e.g., 



2 .2 / in-l)a 



sup r f < 

a<r<b b — a 

Remark 7.5 It follows from the equations above (compare fH] Proposi- 
tion 2.2]) that if there exists r2 > such that 

' p^f{p,x^)dp>{n-l)r2, (7.40) 



then the expansion r(r, x^) will become negative at some value of r strictly 
smaller than r2. If this happens for all x^, then one expects existence of an 
outer trapped surface in the associated space-time. (See [T ^ [35 | li5] for recent 
important results concerning formation of trapped surfaces.) 
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7.4 Determination of Qj^^ 
Recall that we have set 



9ab 1ab = V \ ^Q^^^ J ^AB . (7-41) 



and that (p satisfies the equation 



9ilogy, = -^ = --^, (7.42) 
n — 1 r — 1 

with the initial condition V5(0) = 0. Its integration gives 

ip{r) = r exp ^^dp^ ■ i^AS) 

We assume that the free data Tn and are such that r exists and 
satisfies the conclusions of Theorem 17. 3[ with some a G (0,oo). We have 
then 

1. For small r, using the inequality fl7.2ip . valid for r < a < A^^, 

0<^ <{n- 

we conclude that in such interval we have 
and therefore 

< r - ^(r) <r (^1- exp (^-^^^^^^ ^ < ^^^^^ • (7-44) 

2. For r > a, let ip be as in f l7.19p . we use 

V; ^ ^ < ^ ^ (7.45) 

^2 ^2 ^|]^ _|_ _ ]^^^(.^ j 

to obtain 

^1 log V? = > -— . (7.46) 

r n — 1 l + (ri — IjrCa 
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This shows that if is an increasing function bounded below by 

3. In the case where one assumes that = for r > 6 it holds exactly 

^1 log(^ = > , 4 = (n - l)r-^ - 6 > . (7.48) 
r + df, 

Therefore, using the notation ^p^ := (p{b), 

V{r) = > ¥^6, if r > 6 . (7.49) 

In conclusion if is admissible, and Tu is known and continuous we 
can solve (17. 9p for r on some maximal (possibly angle-dependent) interval 
of r's. Subsequently, (I7.4p can be solved with initial value ip{0) = 0. This 
provides g^B- The quantity ip, hence also Q, depend only on the conformal 
class of 7; the same is true of g, defined by f l7.4ip . 

7.5 Determination of uq 

Once is known we can integrate equation (17. lip for z/q, with the initial 
condition i^o|r=o = 1, 

dii^o = ^(W^i + 'r)uo + i^i^o ■ 

(Note that at this stage any "wave-gauge source" Wi of the form 

Wi{T,gAB,uo,r,x^) (7.50) 

with an appropriate behaviour near r = could be used, though as said 
before, in this section we assume k = 0, and a Minkowski target). The 
equation for z/q reads 

\-g^'rs^B + -}, (7.51) 



i.e., since = — , 



diu' = -^rz/O + ^g^^'rsAB ■ (7.52) 
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This is a linear equation for u^, with coefficients singular for r = 0, but 
continuous for r > 0. Its solution taking a given initial value for tq > 
exists, is and unique for r > ro as long as r and exist and are 
continuous. Note, however, that z/" could go to zero in finite affine time, 
which would lead to a (geometric or coordinate) singularity. 

7.5.1 NCT case 

To study solutions with initial data at r = 0, we start with the 
We have then r = = and fl7.52p reduces to 

d,u' = -^{u' - 1) . 
2r 

The general solution is, for some constant k, 

r\ n — 1 

u — 1 = kr 2 . 
The solution tending to one as r tends to zero corresponds to k 

7.5.2 General case 

To construct a solution tending to 1 as r tends to zero we set 

F := 1 - z/° . (7.55) 
The equation (17.521) for i/° becomes the linear non homogeneous equation 

Y' = -Ky + F , (7.56) 

with F a continuous function (recall the notations r = Trj — ip, = = 
rj^^rsAB and the assumed boundary conditions (14.441) and (I4.49P ) 

F--=lir~ g^^TSAB) ^ ^{(^^^ - g^^'ysAB - ^} = 0(r) , (7.57) 
where V' ^ and 

(r^^^ - g^^ysAB = ^ - m-'^^^SAB. (7.58) 

r 
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NCT case. 
(7.53) 

(7.54) 
= 0, and is 



Incidentally, this implies that F < for initial data such that 

r^fi-V^'sAB > n - 1 . (7.59) 

Now, in the notation of (17. 3p . this can be rewritten in the form 

r^^^-V'^.AB = rV-^[(det7s)^/("-S^T [(det , (7.60) 

such that the two expressions in square brackets have unit determinants. 
Using (p < r, hence r'^ip''^ > 1, the last equation allows one to deduce (17.591) 
from a condition involving only the conformal metric Jab- 

We want to find a solution Y which tends to zero with r, but this solution 
will lead to data for a Lorentzian metric only if remains bounded and non 
zero; that is, if y < 1. 

The homogeneous equation associated to (I7.56P is 

n — 1 1 

Y' = {-^ + -^)Y . (7.61) 

Setting Y = exp Z, this equation reads 

Z' = -!i^4*. (7.62) 

The general solution of (I7.56P is of the form 

Y = wexpZ, with w' = exp{—Z)F . (7.63) 
1. Case < r < a. 

Without loss of generality we can choose 

n-1 , 1 



2 

hence 



logr + ^^ ij{p)dp, (7.64) 



Y{r) = wr 2 exp ^- y ip{p)dpj, with w' = r ^ Qxp (^—- J ip{p)dp^F{r). 

° ° (7.65) 

We find the solution Y := 1 — tending to zero with r (compare (I4.42p ) by 
integrating w' between and r, it gives 



F(r) = r "2' exp (^iy ^{p)dp^ j p"2'exp(^-iy ^lj{x)dx^F{p)dp . 

(7.66) 
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Keeping in mind that there exist numbers C^, a=l,2, such that 
0<ip<Cir , |F(r)| < , 
we see that for < r < a we have 

exp / i^{p)dp^ < exp (^<^i^^), 



^ p'^exp(^-i^ ^{x)dx)\F{p)\dp<C2 dp = ^^^^C2r'^ 

leading to the bound, still for r < a, 



Y < 



n + 3 



exp (^^C'lr^ j 



Since F := 1 — z/° the function z/° is bounded. From f lA.7p the metric will 
have Lorentzian signature, g^^^ being Riemannian, if and only if remains 
bounded and non zero (hence positive since equal to 1 for r = 0). This will 
hold if y < 1, which will be true for any C2 if a is small enough. In vacuum 
C2 is determined by |cTp, so for any a it will hold that Y < 1 for r G [0, a] if 
a is small enough. 

Note that if F < 0, then Y < 0, hence z/° > 1 without restriction on the 
size of a or |cr|. 

2. a < r < 00. 

By the same reasoning as for r < a, the solution Y taking the value Y{a) 
for r = a is 

Y{r) = F(a)+r" V exp (^i ip{p)dp^ J p"^ exp (^-^ ip{x)dx^ F{p)dp . 

3. Suppose that for r > 6 we have g^^ = Lp~'^s^^ ^ hence 

F{t) ^ i{(.^^(r-2 - ^-ysAB - ^} = ^{(1 - rV"')^ - '0} ■ 

We have seen that < r, that is r'^(p~'^ > 1, hence F{r) < and z/°(r) > 
z/0(6). 
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7.6 Vanishing of Hi 

Consider a solution of the wave-map-reduced Einstein equations Rn = Tn 
with initial data on Cq, and with Minkowski target. Suppose that the data 
there satisfy the constraint Ci = Tu. The identity (see (EE])) 



shows that satisfies a linear homogeneous differential equation on Co, 
namely, 

Djf + ^ . (7.67) 

Keeping in mind that D is the covariant derivative of the Minkowski metric, 
in our adapted coordinate system we have 



For all solutions which satisfy uniform bounds near the vertex in 
the {y^) coordinate system, the y'^-components of the wave-gauge vector are 

bounded near the vertex. It follows that H is bounded near the vertex. But 
every solution of (I7.67P which is not identically zero behaves, for small r, as 
j--("-i)/2 along some generators. So, in the uniformly case, we can deduce 
from (17:671) that 



H = , hence also Hi = uqH = . 

Remark 7.6 If we add constraint damping terms as in f l2.13p . we obtain 
instead 

A= (^-^T + em^Hi. (7.68) 

No term proportional to or appears, and hence the damping term is 
compatible with this first step of the wave-map-gauge constraint hierarchy. 
The new term does not change the terms which are singular in r in fl7.67p . 
and hence Hi = is still the only solution with the required behavior. 
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7.7 Scalar and Maxwell fields 



We wish to check that scalar fields lead to equations compatible with the 
required hierarchical structure of the equations. For this, consider a scalar 
field (f) coupled with the gravitational field through an energy-momentum 
tensor of the form 

T^. = d,<pd,<P - + ^(0)) 9o.p . 

In the adapted coordinate system the components relevant for our argument 
are 

Tn = (^10)2^ (7.69) 
Tai = dAjd{4>, (7.70) 

Toi = -uJ:^g^''dA4>dB4> + i^''dA4>di4>-yJ:^t\d{4>? 

-V{ct>)vo . (7.71) 

Keeping in mind that the initial data for the scalar field on Co are provided 
by := 0|co! we see that prescribing provides a Tn which can be used in 
fl73D or in frnHjl (compare ([721) )• 

Next, the relevant components of the stress-energy tensor for the Maxwell 
field Ffj_^ are: 

Tn = V^^'FaiFbi , (7.72) 
Tai = - ly'FAiFoi - V^'^'FabFci + u'' u^^F aiF bi , (7.73) 
Toi = -^uog''''g''''FABFcD-9''''^''FABFci + l^'i^''^''F 

- ^^o^^V'T^i^Bi - ^^°(Foi)2 . (7.74) 

We defer a complete discussion of the Cauchy problem for the Einstein- 
Maxwell equations to separate work. Here we note that if Fia is given on 
the null cone, then fl7.72p is not of the right form for viewing fl7.10p as a 
first order equation for r: Instead fl7.10p should be considered as a second 
order equation for (p, using ( \7A\\ . On the other hand, ( 17.72P is of the form 
(17. 7p . needed for the analysis of the problem when uq has been given. The 
remainder of our analysis of the Ci constraint goes through as before. 
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For further reference, we note that the combination of stress-energy com- 
ponents appearing in the final constraint Cq is 

(7.75) 

for the Maxwell field, and 

f^T^^ + 2^^^T^i + = -g'^'^dA^ (7.76) 

for the scalar field. 

8 The Ca constraint 

The Ca wave-map-gauge constraint operator will be obtained from an analysis 

of _ _ _i _2 

SiA = RlA = RiA + RlA ; (8-1) 

where we have again separated terms including derivatives of Christoffels in 
from the rest in Trivial simplification gives 

Ril = + ^^rf^ - aiF^o - diT^AB ■ (8.2) 

We have by the choice of coordinates F*,*^ = Ff^ = T\i = 0, and therefore 
R-il = ri^(F;^o + Tib) + ^lAi^Bo + ^bi + ^bc) 

~ W'- AO ~ ^ 10^ AB ~ ^ IB^ Al~ ^ IC^ AB- K^-"^) 

We find, for the terms in R^ia, 

doTf; = dog'^iO, lA] + dog'^[B, lA] + ^ly'd^d^, (8.4) 
= ^dAii'^dogu) - i^^XA^dogiB 

+ lii^^id^ - d^UA + 2uBXA'')dmi , (8.5) 

-^iF^o = - ^diiu^idAi^o + dogiA - OiUa)}, (8.6) 

dsTfA - diTAB = dBXA"" - d^iiy'i^BXA'') - dAT. (8.7) 
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And for the terms in R^^ we find 

rLrfs = il'^'idiUA + dAiyo - d^A) - y'^^BXA^'V , (8.8) 

^UKo + + Tbc) ^ XA^'ii^'dsUo + T%c) , (8.9) 
- ^m^M = li'^yidi'yA - dAiyo - d^)d^i , (8.10) 

= ^(z/°)'(9iZ/A + OaUo - d^A - 2iyBXA'')dmi , (8.11) 

- rfor^B = ^i^^XA^idogiB + dlUB - Obi^o - y^ysdogn) , (8.12) 

- TIbTai - Tfc^AB = li^'xA^'id^ - dii^B - Obi^o) - Xb'^^ac ■ (8.13) 

All terms in these formulae can be computed on Cq, except for those that 
contain d^giB or doQu, and whose sum simplifies to 

RiA,ao = -\^^iu'^^) - ir z/°a^ + ^9^(i^°^) . (8.14) 

(We see that all terms quadratic in do derivatives cancel out.) The rest is 
given by 

^lA - RiA,do = ^z^°<9i(9iZ/A - Oa^o - 2ubXa^) + ^bXa^ - uodA{ru°) 

+ ^(9iZ/° + rz/°)(9iz/^ - dAi^o - ^ubXa"")- (8.15) 

8.1 Use of harmonicity functions 

From the identities f l4.27p and f l4.29p we get 

:= 9ab^^ = - ^Ar° + ly^idoQiA + di^A - Oai^o) - 2z/°z/bXa^ + ^a- 

(8.16) 

Hence 

dogiA = -diVA + dA^Q + ^vbXa^ + ^J^i + ^o^a - f a). (8.17) 
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On the other hand 

Ta = Ha + Wa, with HA:=gABH'', (8.18) 
similarly for and Wa- Therefore we have 

dogiA = -diUA + Oa^o + 2i^bXa^ + ^oIa + i^oHa + i'aHi , (8.19) 

with 

fj,:=y%j,W, + WA-fA, (8.20) 
For a Minkowski target, using fl4.18p and f l4.19p . this is 

/a ^ - [x^f^'^scn + ^) + - f ^d) • (8.21) 

8.2 Computation of Ca and Ca 

We see from the identities obtained that -Ri^ is the sum of a linear homoge- 
neous operator Ca on H a '■= VabH and a second order linear operator Ca 
on i>A, both with coefficients depending only on the x^-dependent metric g 
and scalar uq previously determined. (Strictly speaking, ua also appears in 
Ca, but multiplied by H which, with appropriate boundary conditions and 
a Minkowski target, can be shown to be zero at this stage of the argument, 
as explained above.) 

Ca = -\d^(HA + i^Alf) - \r{HA + ^Alf ) + \dA{^jf) . (8.22) 
From dHUD, 

RiA,do = -^i^^didogiA - ^(diu^ + Tv^)dmA + ^5A(z^°5ofi'ii) , 
and using the formula (18.1 9p we find that 

RiA,ao - I^A = —I'^dii-diUA + Oai^o + ^ubXa^ + '^oJa} 
_l(aiz/° + rz/°){-9iZ/^ + dAi^o + ^ubXa"" + ^oIa} 

+ ldA{W, + T). (8.23) 
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Finally, assembling results f lS.lSp and f l8.23p gives 

Ca = RiA -Ca= ^u''di{2diUA - ^r^BXA"" ' ^oIa} 
+ ^{d,u' + ru^){2d,UA - AubXa"" - ^o/a} 
+ ^bXa"" - ^Mr - Wi + 2z/°9iz/o). (8.24) 
It turns out that there is a simple way of rewriting f l8.24p in terms of 

(a ■= -2iy'diUA + ^^''^bXa'' + Y^^^aWi + Wa-Ta , (8.25) 

where f^i := gAud^^^BC (compare fl4.33p ). The vector equals — 2r]^^ 
after using the harmonicity conditions F = W. Note that C,a vanishes when 
va = ^ and = qab- The wave-map-gauge constraint operator Ca can be 
expressed in terms of as 

Ca = -\{diiA + -t^a) + VbXa^ - \dA{r -Wi + 2z/09iz/o) • (8.26) 

Separating different orders of c^i derivatives we get, for a Minkowski tar- 
get, 

Ca = v^'dl^VA - 2u\/drUB + ^° + ^ - ^ W^i - ^"^i^o) d,UA 

- dA Qr + u'd,u, - ]^W^ - \{d, + r){g^Bf''S^^ - f a) . (8.27) 

In the general case, in addition to fl7.6p - (l7.7p one can assume that 

Wa = WA{lAB.^.^a,^A.r,x^) , (8.28) 
TiA = Ti^(source data, 7^5, daAB^ diip, z/q, <9iZ/o, z/^, <9iZ/a, ^ofi-n, r, a;^) , (8.29) 

where as before di denotes derivatives tangential to the light-cone. This is 
clearly compatible with the wave- map gauge f l4.19p . and with scalar fields or 
Maxwell fields as sources (compare Section [TT]). 
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9 Solution of the Ca constraint 
9.1 NCT case 

In the vacuum case with Minkowski target and when = we have Xa^ = 
9ab = r^SAB (therefore f^^ = ^^^), z/q = 1 and r = -Wi = It 
has been shown in [11] that the Ca wave- map-gauge constraint reduces then 
to 

Ca ^ R,A -Ca^ dl.UA + ^-^d,VA + 1 = . (9.1) 

This is a Fuchsian type linear equation, with Fuchsian exponents p = 
and 2 — n. Thus the only solution satisfying f l4.20p . i.e. limr^o^r'^UA) = 0, 

3 — n 

is ua = 0- (In fact, the only solution ua = o(r^~) is zero.) 



9.2 General case 

From the identity f l8.26p we then find 

Lemma 9.1 Assuming ( [7. gp - JT^ and l{8.28\) - l[8.29\) . the wave-map-gauge 
constraint operator Ca = Sia — Ca a first order linear ordinary differential 
operator for the field C,a> with k, as in 1^7. 

1 ~ n — 2 

Ca = --{di^A + r^A) + ^bcta'' - -^^^at - Sak , (9.2) 

where ^a defined as ^8.25\) . which particularizes for a Minkowski target as 

U ■= -2u'd,UA + Au^sXa"" + {W" - VA + lABl''''{S^CD " f ^ 



\B ^ 
CD) 



(9.3) 

Anticipating, we note that va will also appear in the last wave-map-gauge 
constraint Cq through only. 

If one assumes that T^a is known (e.g., in vacuum, or for scalar fields, 
compare fl7.70p ). the homogeneous part of the equation Ca = Tia reads 

-lidiCA + rU = 0, 
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and admits as general solution, keeping in mind that r = (n — l)di log<y9, 

U = L'P-^''-'^ , (9.4) 

for some vector field on the sphere The final solution is of the form 
(19 .4^ . with C,A obtained by integrating the following equation for 

fi — 2 

= 2ip^-'{VBcrA'' -Oat - Oak - T^a} , (9.5) 

n — 1 

with the boundary condition = 0, deduced from the finiteness of 

lim^^ = limr""+^^^ 

r— )>0 r— >-0 

(compare (I4.42p - p.45p . and (l4.49| ). The field va is then obtained by inte- 
grating (19. 3p . with the boundary condition i/^ = at r = 0. These equations 
are a first order linear system of ODEs with coefficients singular for r = 0. 

In the NCT case we have = and W = —x^'g^^ sab, and hence 

2diUA + i^A = , (9.6) 

r 



whose general solution is 



= kAT' 



n — 3 
2 



with kA independent of = r. The solution tending to zero with r, and 
compatible with the boundary condition (I4.43P for n > 2, is = 0. In the 
general case, but with Minkowski target, the equation (19. 3p reads 

n — 3 r f —AR n — 1 4-?/^ r? 

2diVA H - 4:^c(^A + {ri^og sab \ -\i^A = 

r r n — I 

i^olABl'^^'iSED - i^Eo) - ^oU ■ (9.7) 

Setting 

_ n — 3 

T^A = kat- 2 

gives for kA, with fc^ which must tend to zero with r, a differential system 
with coefficients continuous and right hand side tending to zero like , 

dikA = ^kcdA^ + XkA + fiA 
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with 

1 r AR n — 1 4:lb ^ 
A := - {rz/o^^^SAB + ^^}' 

2 r n — 1 

Such a system can be solved by iterated integration starting from k^^^ = 0, 

kf = r{2A:g-^VA^ + Xkt''> + /i^}(p)rfp 
Jo 

The convergence and bound Ik^l < Cr^ results from the bounds of a, A 
and /i. In conclusion, in vacuum, and in the wave-map gauge, the solution 
of fl9.5p exists as long as r does, with ua G and |z^a| < Cr'^. 

9.3 Vanishing of Ha 

The general identity fl2.8p gives in our coordinates 



Ca + Ca = RiA = r[a + ^{^oDaH^ + ^aDiH^ + QabDiRB)^ (9.8) 

with t) the covariant derivative of the target metric, which in this subsection 
will be chosen to be the Minkowski metric, and hence D = D. Therefore, 

( H) — — 

if a metric solves R\^ = Tia and Ca = Tia, we will have, taking H = 
(which, for sufficiently regular solutions, and for a Minkowski target, has 
been justified in Section I7.6P on the left hand side, 

Ca = -^{d,HA + rHA} = ^{uoDJP + uaDJP + g^sDjl^), (9.9) 



For Minkowski target these derivatives are, still for H = 



DaH^ = T%H ^ -rsAsH , D,H^ ^ , (9.10) 



IL ^ Tj^ , T-B -TjC V -BCc^ -Tj , Tj -CdAxB 



D,H^ ^ d,H- + T^cH ^ r'-d.Hc + Hcf'i-^^D ' 2Xz,^). (9.11) 
Therefore (19.91) reads 



2diHA + tHa = {y^rsABt^ - -^A + '^Xa^]Hc • (9.12) 

r 
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Taking leading orders in r near the vertex, as given in Section I4.5[ we 
find 

d,HA + (^S^ + O^)hb = 0, (9.13) 



2r 

where the are 0(r) functions. Hence 



'A 

n — 3 



HA = r-—kA, with ^ifc^ + O^fc^ = . (9.14) 

Standard ODE arguments show that Ha = is the only solution of f l9.12jl 
such that Ha = 0{r), which is the case for metrics having uniform 
estimates at the vertex. 

Remark 9.2 If we add constraint damping terms as in fl2.13p we obtain 
instead, using again H = 0, 

Ca = -\dJiA + \{-r + en^)HA ■ (9.15) 

No term proportional to H^ appears, and hence the constraint damping term 
is compatible with this second step of the constraint hierarchy. The new term 
does not change the leading orders in r of equation f l9.13p and hence Ha = 
is still the only regular solution. 



10 The Cq constraint 

We compute Sqi = Soa^" on Co- We have 

Soi:=Roi-^9oiR, (10.1) 

hence in our coordinates 

= —p^q^'^Rab + RiA^^ - \^^f^Rii ■ (10.2) 

We write 

Rab '■= Rab + Rab 5 (10.3) 



with 



-^AB ^cX'aB ~ Bai ^^AB '■— ^ AB^ a\ ^ AjS^ Ba' (10-4) 
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We will see that the Cq wave- map-gauge constraint is obtained, like the 
other constraints, by decomposing the term in Sqi which has not already been 
computed, 'q^^Rabi into terms defined by data of the degenerate metric on 
the cone and terms which vanish when harmonicity conditions are satisfied 
on the cone. 

Equations f l4.27p and fl4.3ip allow us to express the transversal derivative 
'g'^^doQAB in terms of harmonicity functions, 

+ 2{v'^fv^{-d,VA + vbXa"") + A^'^Vbz/a (10.5) 

^ - r' - d,t' + t\-v'd,u, - ir + + z/V^V^i/A. 

(10.6) 

10.1 Computation of ^^^^^^ 
We have 

Q^^'Rab ■■= g^^'idon^ + d,T\B + OcTab - dATlJ (10-7) 

To compute we proceed in a straightforward way, using the values of the 
Christoffel symbols of the first kind and elementary algebraic relations in our 
coordinates on the cone. Equations (1A.36P and (1A.37P of Appendix |X] are 
useful for the calculations that follow. 

We set 

g^''dor% = h + Ih, (10.8) 

with 

h := g^'^'doll^AB] = -hgAB^OQ^Q^ ^ g'''' OaO^ , (10.9) 

Ih ■■= g''''dog'-[a,AB] ^ {uy^^^g^'' {^d^ - OaUb) - d^r + d^Tc, 

(10.10) 

and 

d^= - i^yO^i - y't'd^i + (^°) V^^, (10.11) 
= {^'f^^d^i - u't^'d^A. (10.12) 



55 



Grouping terms gives 



+ {i^ydogoiT - {v'fv^dog,AT . (10.13) 

We now separate 

g^'^'W^B = Ilh + IVu (10.14) 

with 



II h = 5^^{z/°9i[0, AB] + t^di[l, AB] + t^di[C, AB]} 

(10.15) 



= g^^'d.u^dAUB - ^dogAB) - rd,t' - g^''d^{iy%^)[C, AB]. (10.16) 
Grouping terms gives 

r^aj^^ - ^9i(z/V^a^) -i.V^a^+r^ai(i.°VAi^5) 

-tY^'diXAB-rd.g''- (10.17) 
Finally we have 

r^'dcT^j.B ^ r^'^'^^'dcXAB + rdciv^v') + g^^dct%, (10.18) 
-^^^^aT^, = -^^^9l^(logv^) = -^^^9ls{log(i.oA/ditI)} 

= - g''^{dAi^''dBV^ + i^°9^Bi^o + dAt%c}- (10-19) 

10.2 Computation of Q^^RaI 

We set 

g^^^i^ABK^ - H/^rL} ^ (^2 + ih + iih + iv^ + v, + vh), 

(10.20) 
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with 

h := ^^^r^r^,, ih ■■= ^^^ri#f,. (10.21) 

We find by straightforward computation 

Ih = {^/°VAZ/^-^A^^^^-^'V}(z/°9iz/o + r). (10.23) 

Then we have, recalhng that F denotes Christoffel symbols of the metric g, 

Ilh := ^^^rS^rJ^ ^ (z/V^r + g^^f%){V^j, + u'dci^,) , (10.24) 
Next, 

IV, := -g^^{T\,Tl, + f^if^ J . (10.25) 

Furthermore, 

= -r;^ - v^bXa'' + ^°5a//o, (10.26) 
with {Ca is sometimes called the torsion form) 

tJa = \^^^A^ Ca '■= dogiA + Oai^o - OiUa- (10.27) 

Hence, using f l8.19p . 

G = 29az/o - 2diUA + 2z/bXa'' + MfA + Ha + i^aH"). (10.28) 
In terms of this object we have 

(10.29) 

We set 

v,:= -2^-^^(r°^rgo + ri^rg,) 

= X^''{-('^°) VaCb + 2z/°^^ - 2z/0Vaz/b + 2^1 Vab}- (10.30) 

Finally 

vh - - ^^^rS^r^c = -^^^(-V^xa., + f Sz,)(-°-^xbc + f fc) 

= - (z.°) V^.^Xc'^Xaz. - ^nn-'-^'Xc^ - -9^''f%f^BC- (10.31) 
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10.3 Final computation of q^^Rab 

Adding the results of sections 110.11 and 110. 2[ we get the final result 

+ 2(9i + + r)r' + R- 2g^''v\^T\B - 2g^''V aT\b (10.32) 



R - 2g^^T\^T\B - 2g^^VAT\B (10.33) 



an explicit expression for v\i can be found in f l4.22p ). where 



- 2V^^ = u^dogiA - i^'diUA + 2iy%BXA ' ^^av^ ■ (10.34) 

In this way we have isolated the transversal derivatives in and the 

vector r^^. Decomposing Y a = W a + Ha we find the relation 

- 2T\a = U + Ha + i^'i^aHi, (10.35) 

with defined in ([H25D- 

We note the interesting fact that both the second and third constraints 
naturally break into two first-order equations, with the intermediate variable 
being a Christoffel, respectively F^^^ and g^^Ty^^. 



Remark 10.1 The expression in square brackets in (110.331) can be rewritten 

as 

id,+Tl, + ^)t'+t = -^t'+iy'r^'i^BiyA-l'd^) = g^^T\^. (10.36) 

This shows that 'g^^RAB originally contains only a first-order radial deriva- 
tive of g^^, if we keep the radial derivative of 'g^^d^gAB- It is precisely the 
elimination of the latter object using the harmonicity condition (110. 6p that 
introduces the second-order radial derivative of 'g^^ . Hence, our Co constraint 
operator below will contain such second-order derivative. 

This leads to the following lemmata: 

Lemma 10.2 All terms in'g^^ Rab involving the derivatives dogoi andx^^dogAB 
cancel out. The only new remaining transversal derivative is g^^dogAs, 
which can be eliminated by using T = H +W . 
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Lemma 10.3 It holds that 

Soi = -^uog'^'^RAB + RiAi^"^ - ^i^ot'Ru = Co + Co, (10.37) 

where Co depends only on the quadratic form g on the cone and the W"' while 
jCq is obtained by replacing F by H . 

The explicit formula for Cq reads 

-VaH^ + UH"" + i^ot^djf - u^dAlf 
+ [gooW^ + ^oW^ + uaW^ + ^lyorf' 

-r\[g,^{Tff + 2vjfH' + 2iyATfH^ + 9abH^h''] .(10.38) 

Note that the last line in the previous equation is quadratic in the wave-gauge 
vector H, and equals 



1 ^r/*Tf^ 



^9,.H^H . (10.39) 
All other terms are linear in H. 

10.4 Constraint 

To write the wave- map-gauge constraint Co — Tqi = as an equation for g^^, 
we use the other constraints, which have been satisfied since Ci — La = = 
ii\ = Ha: _ _ _ _ 

R^A = TiA, Ru = Tn. (10.40) 

We find 

-2u%Co - Toi) = 2{di + Kff + 3t(9i + K)f' + {d^r + T^)t^ 

+ 2(^1 + + t)w' + R- ^^^^CaCb + g^^^A^B 

+ t^Tr^ + 2f''T^A + 2f^Toi 
= 0, (10.41) 
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where is the vector fl8.25p and recall that 

«: = z/05iZ/o-^(TFi + r). 

To avoid ambiguities, we emphasise that the right-hand-side of fll0.4ip van- 
ishes identically in wave-map gauge when T^,^ there is replaced by the Einstein 
tensor S^y. This fact reflects the identity, valid for any dimension with our 
choice of coordinates, 

g'^'RAB + ^'':^n + ^f'-S^A + '^f^^oi = 0. (10.42) 
A slightly simplified form of the differential part of the constraint is, using 

mm . 

-2z/0(Co-Toi) = 2(9i + K + r) 

= 0. (10.43) 

Suppose that in addition to (ESD-iO) and (E2H])-(E29]) it holds that 

Wo = Woi7AB,_^,i'A,t\r,x^), (10.44) 
Tn = Tu{...,dogiA,f\dit') , (10.45) 

where ... in fll0.45p denotes the collection of fields already occurring in 08.29 p . 
This is clearly compatible with the wave- map gauge fl4.18p . and with scalar 
fields or Maxwell fields as sources (compare Section [7I7|) . Then (110. 4ip be- 
comes a second order ODE for 'g^^, linear when the vacuum Einstein equations 
and the wave-map gauge have been assumed. 



11 Solution of the Cq constraint 

Throughout this section we assume that the target metric is Minkowski, 
K = 0, and that the relevant components of the tensor T are known (e.g., 
zero). Using the Ci constraint, 

T,, = -(d,T + ^T' + \a\') , (11.1) 
V n — 1 / 
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we find that the Co wave-map-gauge constraint operator can be written as 

(11.2) 

Hence, setting 'g^^ = 1 — a and using using previous notations, the equation 
for the Cq wave-map-gauge constraint, Cq—Tqi = 0, reads as the linear second 
order ODE for a 

L{a) = a" + Ira' + IC^r^ - \a\^)a = $ , (11.3) 
2 2 \n — 1 / 

with $ the function known from the previous sections 
2\n — 1 / 

+Ir- y^'u^B + y^'^A^B - ^'T,^^^ + . (11.4) 

L{a) simphfies to 

r/ \ // 3.n-l ,s , l/(n-l)(n-2) n-2,2 , ,i\ 

L(a) ^ a"+-{——i;)a'+-[^ ^_2^_^V'+-_-V,2-|ar), 

(11.5) 

This hnear equation has smooth coefficients for r > 0, it has a global solution 
with initial data given for r = a > 0. 

We proceed to the study of solutions starting from r — 0. 

11.1 NCT case 

In the NCT case it holds that z/" = 1, = 0, z/^ = 0, r = ^ 

diW^ + tW^ = -{n - l)(n - 2)/r^ = -R and T^p = 0. Hence $ 
Co wave-map-gauge constraint for a = 1 — g^^ reduces to 

3(n - 1) , (n-l)(n-2) 
za -\ a -\ a — (J , 

it is a Fuchsian equation with characteristic polynomial 

2p{p - 1) + 3(n - l)p + (n - l)(n - 2) . 



= 0. The 



(11.6) 
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The zeroes of this polynomial are 

1 — n 

P+ = ' p- = 2-n; 

both negative or zero for n > 2. The general solution of f lll.6p is 

with a± independent of r. The only member of this general solution where 
a tends to zero as r tends to zero is a = 0. 



11.2 General case 

We look for a solution starting from r = and such that 

lim a = lim(r9ia) = . 

r^O r— >0 

We set dia = a' and decompose L as follows 

L{a) = Lo{a) + Li{a) , (11.7) 
where Lq is the Fuchsian operator appearing in the NCT case, 

Lo(a) = a" + —01 + -^a := a" + ^ -oi + ^ -a , (11.8) 

r r 2r 2r 



3 r {n-2) \ (n-2 



a\ \ > a. 



and 

Li(a) = a\a + h\_a := -T^V^a' + \ -— — —'il) + ^ ( - — '^i)^ 

2 I r 2 V 72 — 1 

In order to use the idea of the Fuchs theoreno, we write the second order 
equation flll.3p as a first order system for a pair of unknowns v := (fi,f2) 
by setting v\ := a, V2 = ra', hence ra" = — r^^t>2 + The system reads 

rv[-V2 = , 

rv'2 + (flo - 1)^2 + &o^i + r{aiV2 + 6iri;i - r$) = . 
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See e.g. [9], Appendix V. 
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It is of the form 

rv' + Av^r{Fi{r)v + Fo{r)}, (11.10) 
with A the constant hnear operator 

^-(! ~\) 

\bo oo - 1 / 
whose eigenvalues are found by solving the equation 

det ( ° ~ ^ ~} ) = + 1^(1 - ao) + 6o = . 

The solutions are the opposites, — p±, of the characteristic indices computed 
in the NCT case, hence nonnegative. Further 

Fi{r)v= ( ° , 1, Fo{r)= f M , 

where 

are bounded functions smooth away from r = 0, as well as r$. What has 
been said shows that Fq and Fi are continuous at r = for admissible 7^^- 

Lemma 11.1 Let 

rv' + Av ^ r{Fi{r)v + Fo{r)} , (H-H) 

be a linear differential system with A a constant linear operator with non- 
negative eigenvalues. Let Fi be a continuous linear map and F2 a continuous 
function, for < r < Tq. The system admits one and only one solution in 
C^{[0,ro\) which vanishes at r — 0. 

Proof. We set v = Mw, with M a 2 x 2 matrix satisfying the homoge- 
neous equation 

rM' + AM = . 

We choose for M the matrix 
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The equation flll.lUp reads then 

w' = r^{Fi{r)v + Fo{r)} . (11-12) 

Hence the equation (111.111) together with the condition v\r=o = are equiv- 
alent to the integral equation 

vir) = f\T-'p)''{F,{p)v{p) + Fo{p)}dp . (11.13) 
Jo 

We have 

sup \{r-^p)^\ < 1 . (11-14) 

0<p<r 

We set, with a an arbitrary positive number, 

Ci ■= sup |Fi(r)| , Co := sup |Fo(r)| . (11.15) 

0<r<a 0<r<a 

The integral equation fill. 131) can then be solved by iteration, setting 

<r)= ! {r~^p)^F^{p)dp. 
Jo 

Hence for r < a 

\vo{r)\ < rCo , 

pr 

Vi{r) := / {r-^p)'^Fi{p)vo{p)dp + Vo{r) , 
Jo 

and so 

\vi{r)-Voir)\<^r^C,Co, 

pr 

Vn+i{r) := / {r'^p)'^Fi{p)vn{p)dp + Vo{r) , 
Jo 

pr 

\Vn+l{r) - Vn{r)\ < / Ci\Vn{p) - Vn-l{p)\dp . 

Jo 

Assume that, as satisfied for n = 1, 

Then the same inequality is satisfied when replacing n by n+1. The sequence 
Vn converges therefore in uniform norm to a limit v, solution of the integral 
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equation f lll.l3p :. hence of the differential equation f ill. lip . This solution 
V := (a,rQ;'), a(0) = (ra')(0) = 0, is defined, continuous and bounded for 
any finite r. ■ 

We deduce from this lemma the following theorem. 

Theorem 11.2 In the interval of r > 0, possibly angle dependent, where 
the Ci constraint has a global solution and pq > 0, the Co wave-map-gauge 
constraint with coefficients deduced from the solutions of the Ci and Ca con- 
straints admits a solution g^^ = 1 — a with a(0) = (r«')(0) = 0, a and ra' 
which are in r. The solution is global when it is so of the solution of the 
previous constraints, since the system is linear. 



11.3 Vanishing of Hq 

In previous sections we have seen how to achieve H = H =0, and hence 
Ci = Ca = 0. Specializing equation (110. 38p to this case we get (with Hq : = 

Co = -djio + ^(Wi - t)Ho . (11.16) 
On the other hand, the identity 

S'oi = 4f + \{9^aDiH^ + gic.DoH'' - goiDaH'^) (11.17) 
reduces on Cq to 

Soi = + ^{gooDjp + uaDiH^ - uoDaH^)- (11-18) 

A 

Using again the conditions H = H = we have, for an arbitrary target 
space in adapted coordinates, 

biH^ = , DiH"^ = , DaH"^ = f^i# , (11.19) 
and hence flll.lSp further reduces to 

Soi^siT -IriJlo. (11.20) 

For a solution of the Einstein equations in wave-map gauge it holds then 
that ^ 

•S'oi = Co + Co = Toi — -T'^iHo- (11-21) 
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Therefore, when H = H =0 and the initial data satisfy the wave-map- 
gauge constraint 

Co-Toi = 0, (11.22) 

then Ho satisfies the equation 

^1:^0 = ^(Wi - r + ri,)Ho . (11.23) 
For a Minkowski target, and using the boundary conditions ( I4.20p . we have 
f^^ 1 , \im{rWi) = hm(-rV^SAB) = -(^^-1), hm(rr) = n-1 , 

r r— >-0 r— >0 r— >-0 

(11.24) 

and hence equation flll.23p takes a Fuchsian form 

rdiHo + ^^^^0 + rMHo = . (11.25) 

with M a continuous function up to r = 0. 

We want to prove that Hq = when the spacetime metric is a solution 
of the Einstein equations in Minkowski- wave- ap gauge; in this case the wave 
gauge vector H is C^, then Hq tends to a finite limit at the vertex. The 
equation flll.25p implies that this limit is zero hence that the only solution 
is zero. 

Remark 11.3 If we add constraint damping terms as in (12.131) we obtain 
instead, using again H = H =0, 

Co = -djio + ^(Wi - r + epn^)Ho . (11.26) 

This new term does not change the leading orders in r of equation f lll.25p 
and hence Hq = is still the only regular solution. We conclude that the 
addition of constraint damping terms is fully compatible with the wave-map- 
gauge constraint hierarchy. 
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12 Wave-map gauge constraints: a summary 



We have defined Co to be the cone represented in R""*'^ by the Minkowskian 
cone 

2/° = r, r':=Y,{yl\ (12.1) 

i=l 

equivalently 

a;0 = 0, x^=r, Q\x^) = -. (12.2) 

r 

We have considered on Co a non degenerate quadratic form given in 
coordinates by: 

'g^Q^dx^Y + 2vQdx^dx^ + 2vAdx^dx'^ + g^s^^dx^dx^. 

We have proved (recall that admissible means hypotheses on smoothness and 
limits at the vertex spelled out in various sections) 

Theorem 12.1 1. Let g he a given admissible degenerate quadratic form on 

> 

g = gABdx^dx^. 

There exists on Cq, for some T with < T < Tq, coefficients z/q, i'a, goo sat- 
isfying the vacuum Einstein wave-map gauge constraints, and unique modulo 
admissible vertex limits. 

2. An admissible degenerate quadratic form g together with a non vanish- 
ing z/q can be determined on Cq, for some T with < T < Tq, from the first 
vacuum Einstein wave-map gauge constraints, an admissible quadratic form 
7 and a scalar function k being arbitrarily given. Then g is conformal to 7 
and depends only on its conformal class, uq is linked to k by the differential 
equation 1^7. They are unique under admissible vertex limits. 

When g is known, ua and are determined as in point 1. by the second 
and third wave-map gauge constraint and admissible vertex limits. 

13 Local geometric uniqueness for the vac- 
uum Einstein equations 

In this section only the vacuum Einstein equations will be considered. 
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Recall that two spacetimes (Kjfi'a) and are considered as (both 

geometrically and physically) the same if there exists a diffeomorphism : 
K — )■ H such that on Va it holds that Qa = (l>*9b- We have said before 
that given a C^'^ metric Qa on a manifold Va and Oa G \4 there are in 
some neighbourhood of Oa normal coordinates centered at Oa, where the 
characteristic cone Cq^ is represented, for < ?/° < by the equation of a 
Minkowskian cone in 

n 

1=1 

The null rays issued from Oa are represented by the generators of this cone. 
We have defined adapted null coordinates by setting 

n 

xl:=ra-y'a, xl = ra, with rl = Y,{yi^f (13.1) 

1=1 

and x^ local coordinates on the sphere S'"'^^. In the coordinates x° the metric 
Qa reads on the cone Cq^ 

9a,oo{dXaf + 2vafldxldxl + 2va,Ad^ld^a + QaABd^adx^. (13.2) 

We have shown moreover (see Section 14. 5p that it is possible to choose the 
coordinates so that it holds 

hmr-^{ga,AB-rlsAB) =^ , \i-mr-'^di{-g AB-rlsAB) = ^ , (13.3) 



limr„ ^(z/a,o - 1) = limr„ Va,A = limr^ ^^iZ/q^a = , (13.4) 

r— >-0 r— >-0 r— >-0 



and even 



linir^ do{ga,AB - rlsAs) = , (13.5) 



while 



lini r„ ^doga,iA = hm ^dogaflA = , (13.6) 

1 — ^0 r— >-0 



lirnr^ ^(^a,oo + 1) = , lim^i^ qo = = \im doga,oo ■ (13.7) 



Having chosen such coordinates y^, respectively y^^, for the metrics ga and 
(7b, we obtain a diffeomorphism 0^ by y^iy^) '■= Ua^ defined in the subset 
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l/i^T := min(Ta, Tfe). Such a diffeomorphism will be called canonical. We 
remark that canonical diffeomorphisms are not unique, and that above we 
have not required r to be an affine parameter. 

The metrics gh and (t)M,*9h ^^re geometrically equivalent, and one has equal- 
ity of components {(pN *9b)^^{ya) = db^iyb) for = ?/°. The coordinates 
are normal for both metrics and they satisfy in the coordinates the vertex 
limits fll3.3f[TX7|) recalled above. 

To study the geometric uniqueness of our characteristic Cauchy problem 
we first consider two metrics Qa and gh on the same manifold which satisfy 
the characteristic Cauchy problem on the same cone Co- We will prove the 
following theorem, using the notations given in previous sections for Co and 
Yo (note that we are not assuming an affine parameterisation of the cone 
generators here): 

Theorem 13.1 Consider two smooth solutions Qa and gt in Yq of the Cauchy 
problem for the vacuum Einstein equations Ricci{g) = with data on the 
cone Cq, characteristic for both metrics. There exists T' <T so that ga is 
equivalent to gb in Yq if and only if they induce on Cq the same degenerate 
quadratic form satisfying in the coordinates x° the vertex limits U3.3[\T371\). 



Proof. We put the metric ga in Minkowski wave- map gauge by con- 
structing a wave map fa, that is a solution of the semilinear, tensorial, partial 
differential equations which read in abstract index notation 

□ 9.M: = ga'idlX - ra,l, + d^f:d,fPa^%) = , (13.8) 

which on Cq is the trace of the identity mapping / of R"'^^. To simplify the 
writing we suppress the index a in the following computations, valid for any 
metric g with normal coordinates and adapted null coordinates x", we 
will reestablish a and b in the conclusions. 

The components /" and of the image point are linked by the same 
relations as the coordinates y and x. They take in coordinates the initial 
data 

7° = 0, f=x\ T = x^, fora;° = 0; (13.9) 
and in the coordinates the initial data 

T = y\ f = r, (13.10) 
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we see that in the y coordinates the initial data are the trace on Cq of the 
set of C°° functions on R""*"^ 

P = y\ L' = y'. 

The existence of a wave map / in some Yq"" taking these initial data 
can therefore be proved using the Cagnac-Dossa theorem. In fact, since 
the equations are linear in a coordinate system where the F's vanish, the 
usual linear theory p2] suffices to obtain the result. The resulting wave map 
extends to a mappinj^. 

To prove that it is a diffeomorphism at least in a neighbourhood of the 
vertex we first remark that our definitions imply 



dj' = 61 (13.11) 

To study the derivatives do we return to the x coordinates and consider the 
set of functions 

f-x\ f^-x^. (13.12) 

They vanish on Co, so do therefore their tangential derivatives on Co, hence 
by application of the Lemma 14.21 



limao/i = 0, \imdAp = 0, limao/o = l. (13.13) 

i — ^0 r— >-0 r— ^0 

By definition of the coordinates x and y we have 

r = (s(f)2)i-r 

hence 

M!-|o/!--|./°-^/°-«o/° (13.14) 

while depends only on and f^. Therefore 



limaof = 0, lim9o/o = l. (13.15) 

r— >0 r->0 

Since the Jacobian of the C^ mapping / tends to 1 at O, it is a diffeomor- 
phism, between at least small neighbourhoods of O. 



^though not in general to a wave map. 
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The initial data, trace g^^^ of the metric g^^^ in wave gauge are hnked 
with the original g by the classical relation 



9a0 = ^a/'^/^M? • (13.16) 

The values of dif in the coordinates show the equality of quadratic forms 
g^^'^ = g, indeed in these coordinates: 

^Sf ^ = ^11 = , gi^^ =9iA = , gfi = 9ab ■ (13.17) 

Since g^^^ is in wave gauge, and satisfies the vacuum Einstein equations, 
its trace 'g^^^ satisfies the wave- map gauge constraint Ci = 0, and z/q^"* satis- 
fies the same equation than z/q, 

^i^^o = '^o { + T^i^o ^ + ^) + + 7}> (13.18) 

r / T n — i 

since the coefficients depend only on ^, to show that z/q tends to 1 at the 
vertex like z/q we use the identity 



Z.0 = dof'dj'^gif . (13.19) 
and the limits ffm^ - ffm^ give 

1 = lim z/q = lim ui!^"' . 

Uniqueness of solutions of (113. ISp (with non zero r) with this limit at O shows 
that the function z/q^^ depends only on g, i.e. u^^J = z/[^^ since ga = gb, that 

is ga,AB = 9b,AB- 

As a consequence, the wave-map-gauge constraints Ca,A = written with 
g^^ = g^^j^ and u^J^q are the same equation for ly^^J, a = 1 or 2. The vertex 
limit of z/)^ will be deduced from the definition: 

i^A := VoA ^ dj^d^gif = dJ^vTx (13-20) 
which implies using (113. 9p (compare (I4.43P ) 

= limr~^z/^ = lim r'^z/^^'', limz/"^ = lim'q^^i'B = limr'^s^^UB = . 

r^-O r-5.0 r^O r^O r-s>0 

(13.21) 
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Differentiating (113. 20p gives 
We have 

\iYa.r~^diUA = lim (r^'^didof^gf^) + lim {r^^dof^d{gf^] 

r— i>0 r— i>0 r— >-0 

with, by and (fmsj) . 



r— >0 

and 



1 — >0 r— 5-0 



hm(r ^didof^gf;)!) = \im{rdidof^)sAB- 

r— )-0 r— i>0 

Taking the trace on the cone of the wave map equation, with Minkowski 
target, gives 



-2^^'^(rf^-fi^B)-^'^''(rL^ + ric--f^c) = o. (13.22) 



We have 

1. ^ 



hm FiQ = and hm dof^ = , 

r— )-0 r— s>0 

hmrffo = hm ^{-rz/^z/°9ofi'ii + rg^^{dQgiB + 9iZ/b - 83^0)} = . 
FinaUy, for a wave map / 



hm (rdidof^) = 



hence, since Mmr^or ^OiVa = 0, we obtain 

\im Ot u^f^ = . 

r-fO 

Since i^a.A and u^^a satisfy the same equation and the same boundary condi- 
tions, they are equal. 

It remains to analyse the boundary conditions for the functions ^^00; 
which again satisfy the same equation for a = 1 or 2. We have 

900 = dof%Fg[f. 
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It implies: 



r->-0 r— >-0 

The previous limits imply then 

lim^oo = lim^o?^ = -1- 

Also 

hence using previous limits 

linir9i^oo = li™;{r(9i^[,{^) + 2r {did^Pg^^^ + 9i(9o/V[,^))}. 

We have, by definition of a wave map with Minkowskian target, 



Hence 

limr{2(9i9o/o - limr") + ^ } = 

r->o r r 

which gives 
One finds also 



limr9i9o/° = 

r— s>0 



limrc^i^o/^ = 



hence 

limr^i^oo = limr9i^S)J) = • 

r— >-0 r— >-0 

We have proved that % = gt on Cq implies 'g[^^ = 'gl^^ on Cq hence by 
uniqueness for the hyperbolic system of the Einstein equations in wave gauge 
gi^^ = g\^^ in Yq . The metrics ga and gh are geometrically equivalent. 

The reverse implication is trivial. ■ 

Our next result, one of the main results of this paper, is a straightforward 
corollary of Theorem 113.11 
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Theorem 13.2 Given points Oa G Ki (^nd G H denote by Cq^ o.f^d Cp^ 
the characteristic (null) cones of smooth Lorentzian metrics Qa on Va and 
on Vb- Denote by J+ the future of the point Oa in the metric Qa- There are 
neighbourhoods Ua of Oa and Ub of Ob such that the spacetimes {UaH , ga) 
and {Ub n Jb,gb) are locally geometrically the same if and only if the pull 
back 4>*j^gb, where (p^ is a canonical diffeomorphism of Ua onto Ub, equals (ja- 

Proof. The spacetimes iUbnJ+,gb) and (0^^(f/fenJ+) C [/„n J+), 0^c/fe) 
are geometricaUy equivalent. Theorem 113.11 shows that the second one is 
locally geometrically equivalent to (t/aH J+), (7^), the conclusion follows from 

the fact that 4>*j^gb = (t)*N9b and satisfy the required vertex limits. ■ 

From the Uniqueness Theorem 1 1 2 . 1 1 for the constraints one deduces straight- 
forwardly a formulation of geometric local uniqueness starting from data 7 
and K. 

14 Conclusions, and open problems 

We have shown that the trace ^ on a characteristic cone of a solution of 
Einstein equations which is also a solution of the reduced Einstein equations 
in wave-map gauge satisfies necessarily a set of n -|- 1 equations which we 
have called wave-map gauge constraints, written out explicitly and solved. 
We have shown that, conversely a solution of the reduced Einstein equations 
in wave-map gauge with trace satisfying these wave-map gauge constraints 
satisfies the original Einstein equations. Finally we have shown that every 
solution of the vacuum Einstein equations is locally (i.e. in a neighbourhood 
of the vertex) isometric to a solution in wave map gauge, uniquely deter- 
mined, modulo some regularity conditions, by the degenerate quadratic form 
induced on the characteristic cone by the spacetime metric. 
There remain many interesting open problems: 

• Determine the minimum regularity, in particular at the vertex, under 
which the initial data lead to a local solution (see also [12J). 

• Extend our analysis to a characteristic cone with vertex at (cf. [21]). 

• Study the asymptotic behaviour of the solutions of the wave-map gauge 
constraint equations at future null infinity. 
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• Prove global existence for small initial data of solutions of the Einstein 
equations in higher dimensions by a conformal method, as was done 
for the spacelike Cauchy problem with data identically Schwarzschild 
outside of a bounded region \10\ . 

• Prove global existence using the approach of Lindblad-Rodnianski [3H1 
[39] (compare [Il[26]). 



A Collected formulae 

The metric on Co = {x^ = 0} is written as 

g = gQQ{dx'^y + + + 'gj^^dx^dx^ , (A.l) 

and recall that we do not assume that this form of the metric is preserved 
under differentiation in the x°-direction. Here and elsewhere we put overbars 
on the relevant quantities whenever the formulae hold only on Co- The 
inverse is 

g^ = t'df + 2t^d,dA + 2f'dodi + g^'^'dAdn , (A.2) 

with 



(A.3) 



r = ^, = -rg^^B , r = irn-9oo + 9^%a9ob 

We introduce the special notations 

1^0 ■= 9oi , i^A ■= 9oA > 9 ■= VAB^x^dx^ , (A.4) 

^° := r = -, := ~%it^ = g^^'^B . (A.5) 

Then 

-lA ^ _^0^A -11 ^ (^0)2(_- ^ ^ ^A^^^ (A.6) 

The determinant reads 

a/I det^l = z^oVdet^s ■ (A.7) 

The Levi-Civita connection of the metric g^B will be denoted as Va, with 
corresponding Christoffel symbols P^^ with respect to the derivative d^- 
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We have the following Christoffel symbols on the null hypersurface, 



^ 00 




— V^( — (9i f7n^> + 2(9nf7ni 1 

2*^ V "-^lyoo ~ '^"^oyoiy ) 




^ 01 




'-'oyii ) 


(A 9) 


11 
on 




, 

^I/°I/^(9AOno — 29oOoa) + ^^^H— (?1^00 + 2(?oOoi) 

+^i^°5o5'oo , 


(A.IO) 
(A.ll) 


r' 

^ 01 




2 V lyoo ~ ^^i^A^^o "^I'^A •^oyiA^ ~y <^oyiiy ^ 


(A.12) 


i 11 








r° 

^ AO 




— ( f) A Urt -\- f)nni A — (9i Z/ 4 1 
2*^ V A ~ "-'OyiA '^I'^AJ 1 


(A.14) 






, 


(A.15) 


^ AO 




z 

+ ^^^^(^Ai^0 + ^OS'IA - QiVa) , 


(A.16) 






^I/°(aAt^O - ^OS'IA + diUA - l^^diQAB) > 


(A.17) 


^ AB 




-\^^d{gAB , 


(A.18) 


^ AB 




1 ■ ■ 1 

-1/ [VaI^B + ^bT^A - OogAB) - OiQaB ) 




^ 00 




~2^ ^^9oo +2 '^iS'oo + fl' O'ofi'OA — T^T^ Oogol , 




^ 01 




2^ I'^oyiA t^ii^A tJ^Ai^oJ 2 ^oyii ; 




^ AO 




0, 

--V^V^idoQiA + 5^i^0 - ^li^A) 

+ ^^^^(VAi^B - Vfii/A + OgQab) , 


(A.22) 
(A.23) 






^g^^diQAB , 


(A.24) 


^ AB 
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(A.25) 



The remaining ones are obtainable by symmetry. Note that in spite of hav- 



ing Qj^Q = qab, the Christoffel symbols F^^ (a part of F^j^) and F^^ (the 
Christoffel symbols of Qab) do not coincide in general. 
We note the following traces of the Christoffel symbols: 

F^^ = u'd^^ + ^t'd^,-u%^d^+^g^''d^, (A.26) 
% = i^'d,vo + \g^^d{g^^, (A.27) 
r:;, = y'dAyo + lr^'dA-gBc- (A.28) 
The harmonicity vector on the null surface reads: 

r° = {v'fd^i-\v'9^''d{gAB. (A.29) 
F^ = z/V^Vb^^^ + t'y'd,uo - \t'9''''d{g^s + {u'fu^u^d{gAB 

,0n2o- o/',,0n2,,Ao 1, ,0^7^50—; I -'-,,0^111 



H^yd{g,, - 2{uyu''d,UA - -y'r^'d.gAB + -y't'd^g,, (A.30) 

- '9i('^o^'Vdet^s) 1 0-AB7T-— , 1 n_ii^T— 

- y g y B^A / , ^ 7;^ g o^gAB + t;^ g oogn , [A.6L) 

z/ovdetfifs 2 z 

i - -g u ObI^O + g ^ CD+ 2^ ^ ^ OigBC "5' ^ ^ Oig BC 

+z/0(^^^9iz/B + g^^'d^ - z/V^a^) , (A.32) 
^o^r = -v\^dAVo + g^''dBVA + v'^d{g^^+t^d^v^-v%^diVA 

—g^^'d^ + z/°z/^^ - ^t'd^i , (A.33) 

gi,T' = -^^''^9i^AB + z.°a^, (A.34) 
^^/'^ = -u\dAiyo-d,iyA-d^+iy''d,gAB)+g''''gAD^Bc- (A.35) 

Q 

(In the main body of the paper we also use F^ := ^^^F , see f l4.35p .) 
The following formulae are often used in our calculations: 



d,gOo ^ -{uydogu , dog^^ ^ -y\-i^'i^''d,g,,+g''^d,g,c) , (A.36) 



dog'' ^ -{{uydogoi + ly't'dogn - (^")V^9o(7ic)} • (A.37) 
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The scalar wave operator acting on a function / reads 



+9^''dAdBf . (A.38) 

The tensor computations in this article have been checked with the com- 
puter algebra framework xAct \A0\ . 
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